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STABILITY  CRITERIA  FOR  DIFFERENCE  APPROXIMATIONS  TO  HYPERBOLIC 
SYSTEMS,  AND  MULTIPLICATIVITY  OF  MATRIX  AND  OPERATOR 


Principal  Investigator:  Moshe  Goldberg 


ABSTRACT 


Research  completed  under  Grant  AFOSR  -83-0150  by  Moshe 
Goldberg  consists  of  the  following  two  topics: 


9— (afConvenient  stability  criteria  for  difference  approximations  to 
hyperbolic  initial-boundary  value  problems,  -5.) 


£z— fb^uKiplicativity  and  stability  of  matrix  and  operator  norms. 


^  j  °  /T  .  \ 


J 
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STABILITY  CRITERIA  FOR  DIFFERENCE  APPROXIMATIONS  TO  HYPERBOLIC 
SYSTEMS,  AND  MULTIPUCATIVITY  OF  MATRIX  AND  OPERATOR 

Principal  Investigator:  Moshe  Goldberg 


1.  Convenient  Stability  Criteria  for  Difference  Approximations  to  Hyperbolic 
Initial-Boundary  Value  Problems 


Consider  the  first  order  system  of  hyperbolic  partial  differential  equations 


9u(x,t)/9t  *  A9u(x,t)/9x  +  Bu(x,t)  +  f(x,t),  x £ 0,  t £ 0,  (1  .la) 


where  u(x,t)  »  (u(1)(x,t), ....  u(n,(x,t))'  is  the  unknown  vector  (prime  denoting  the 

transpose),  f(x,t)  »  (f<1)(x,t) . ^(x.t))’  is  a  given  n-vector,  and  A  and  B  are 

fixed  n  x  n  matrices  such  that  A  is  diagonal  of  the  form 


A 


0  A" 


a1  o  a«  n 

A  >  0,  A  <  0, 


(1.2) 


with  A1  and  A11  of  orders  k  x  k  and  (n  -  k)  x  (n  -  k),  respectively. 

The  solution  of  (1.1a)  is  uniquely  determined  if  we  prescribe  initial  values 

u(x,0)  -  u(x),  x  £  0,  (1.1b) 

and  boundary  conditions 

u "(0,t)  »  Su'(0,t)  +  g(t),  t*0,  (1.1c) 

where  S  is  a  fixed  (n  -  k)  x  k  coupling  matrix,  g(t)  a  given  (n  -  k)-vector,  and 
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u 


(U 


(1) 


A-. 


u 


<U 


(k+1) 


u(nV 


(1.3) 


a  partition  of  u  into  its  outflow  and  inflow  components,  respectively, 
corresponding  to  the  partition  of  A  in  (1 .2). 

In  the  past  five  years,  E.  Tadmor  and  I  [22-24]  extended  our  previous  results 
in  [19,20]  to  achieve  versatile,  easily  checkable  stability  criteria  for  a  wide  class 
of  finite  difference  approximations  to  the  above  initial-boundary  value  problem. 

More  specifically,  introducing  a  mesh  size  Ax  >  0,  At  >  0,  such  that  X  *  At/Ax 
*  constant,  and  using  the  notation  vy(t)  *  v(vAx.t),  we  approximated  (1.1a)  by  a 

general,  basic  difference  scheme  -  explicit  or  implicit,  dissipative  or  not,  two- 
level  or  multilevel  -  of  the  form 


Q.1vy(t  +  At)  =  ^  Qovv(t  -  oAt)  +  Atbv(t),  v  *  r,  r  + 1 . 

am  0 

Q  »  2j  A.  E*.  Ev  a  V  ,  O  a  -1 . S, 

a  jo  v  v  +  1 

J”-f 


d-4) 


where  the  n  x  n  coefficient  matrices  Aj0  are  polynomials  in  XA  and  AtB,  and  the 
n-vectors  bv(t)  depend  on  f(x,t)  and  its  derivatives. 

The  difference  equations  in  (1.4)  have  a  unique  solution  v  (t  +  At)  if  we 
provide  initial  values 

Vv(pAt)  -  v(pAl),  p  a  o . s,  v  a  0, 1,  2 .  (1.5) 

and  specify,  at  each  time  level  t  a  pAt,  p  a  s,  s  +  1 . boundary  values 

vv(t  +  At),  v  a  o,  ....  r  - 1.  Such  boundary  values  are  determined  by  boundary 
conditions  of  the  form 
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T^VJt  +  At)  -£l*”vy(t  -  oAt)  +  Atdv(t),  v  *  0 . r- 1, 


o*o 


T<v|  -  f  CMe‘, 

o  r“*  )o 

l-o 


(1.6a) 


&  *  *1  f  •••!  Q| 


[yj 

where  the  n  x  n  matrices  C.  depend  on  A,  AtB  and  S;  and  the  n-vectors  d  (t) 

&  v 

are  functions  of  f(x,t),  g(t)  and  their  derivatives. 


Our  intention  was  to  interpret  the  difficult  and  often  stubborn  Gustafsson- 
Kreiss-Sundstrdm  (GKS)  stability  criterion  in  [26]  in  order  to  obtain  simple  and 
convenient  stability  criteria  for  approximation  (1.4)  -  (1.6a).  While  we  were 
unable  to  meet  this  goal  for  general  boundary  conditions  of  type  (1.6a),  we 
managed  to  achieve  rather  satisfactory  results  under  the  further  assumption 

(v) 

that,  in  accordance  with  the  partition  of  A  in  (1 .2),  the  C;  can  be  written  as 


C,II(v) 

jo 


c«  »<v) 
jo 


(1.6b) 


where 


the  c' 1  are  independent  of  v, 
jo 

the  c' 1  are  diagonal  when  B  *  0, 

|o 

the  C1  ll(v>  ■  0  when  B  ■  0, 

JO 

CH  ll(v)  *  0  for  j  >  0  and  a  >  -1  when  B  =  0. 
jo 
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(1.6c) 

(1.6d) 

(1.60) 

(1.6f) 


The  essence  of  (1.6c)-(1.6e)  is  that  for  B  ■  0,  the  outflow  boundary 
conditions  are  translatory  (i.e.,  determined  at  all  boundary  points  by  the  same 
coefficients),  separable  (i.e.,  split  into  independent  scalar  conditions  for  the 
different  outflow  unknowns),  and  independent  of  outflow  values.  Assumption 
(1 .6f)  implies  that  for  B  -  0,  the  inflow  values  at  the  boundary  depend  essentially 
on  the  outflow. 

It  should  be  pointed  out  that  our  outflow  boundary  conditions  are  quite 
general,  despite  the  apparent  restrictions  in  (1 .6c)-(1.6e).  Indeed,  (1.6c)  is  not 
much  of  a  restriction,  since  in  practice  the  outflow  boundary  conditions  are 
translatory.  In  particular,  if  the  numerical  boundary  consists  of  a  single  point, 
then  the  boundary  conditions  are  translatory  by  definition,  so  (1.6c)  holds 
automatically.  The  restrictions  in  (1.6d),  (1.6e)  pose  no  great  difficulties  either, 
since  they  are  satisfied  by  all  reasonable  boundary  conditions,  where  for  B  »  0 

the  usually  reduce  to  polynomials  in  the  diagonal  block  A1,  and  the  cJJl(v) 
vanish. 

We  realize  that  in  view  of  the  restriction  in  (1.6f)  our  inflow  boundary 
conditions  are  not  quite  as  general  as  the  outflow  ones.  They  can,  however,  be 
constructed  to  any  degree  of  accuracy  (see  [20]);  and  if  the  boundary  consists  of 
a  single  point,  then  such  conditions  can  be  achieved  in  a  trivial  manner,  simply 
by  duplicating  the  analytic  condition  (1.1c),  which  gives 

v[[(t  +  At)  -  Sv|j(t  +  At)  +  g(t  +  At). 


Throughout  our  work  we  assume,  of  course,  that  the  basic  scheme  (1 .4)  is 
stable  for  the  pure  Cauchy  problem,  and  that  the  other  assumptions  which 
guarantee  the  validity  of  the  GKS  theory  in  [26],  hold. 

The  first  step  in  our  analysis  was  to  reduce  the  above  stability  question  to 
that  of  a  scalar,  homogeneous  problem.  This  is  obtained  by  considering  the 
outflow  scalar  equation 

9u/9t  ■  adu/0x,  x£0,  t£0,  a  »  constant  >0,  (1.7) 

for  which  the  basic  scheme  (1 .4)  reduces  to  the  homogeneous  scheme 
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Q-1vv(t  + At) 


tQ  v  (t  -  crAt) 

O  V 

o-0 


Q 

o 


a  =  -1, s, 


(1.8a) 


and  the  boundary  conditions  (1 .6)  reduce  to  translatory  conditions  of  the  form 


T.^t  +  At) 


Z,  Vv(t ' oA1) 

o-0 


T 

a 


a  »  -1 


q. 


(1  -8b) 


where  a  and  c.  are  scalar  coefficients. 

jo  JO 

Referring  to  (1 .8)  as  the  basic  approximation,  we  proved: 

Theorem  1.1  [24],  Approximation  (1.4)-(1.6)  is  stable  if  and  only  if  the 
reduced  outflow  scalar  approximation  (1.8)  is  stable  for  every  eigenvalue 
a  >  0  of  A1.  That  is,  approximation  (1.4)-(1.6)  is  stable  if  and  only  if  the 
scalar  outflow  components  of  its  principal  part  are  all  stable. 

This  reduction  theorem  implies  that  from  now  on  we  may  restrict  our  stability 
study  to  the  basic  approximation  (1 .8). 

In  order  to  introduce  our  stability  criteria  for  the  basic  approximation,  we  use 
the  coefficients  of  the  basic  scheme  (1 .8a)  to  define  the  basic  characteristic 
function 


P(z,k) 


± 


-o- 1 


a.  z 


JO 


id. 
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Similarly,  using  the  coefficients  of  the  boundary  conditions  in  (1 .8b)  we  define 
the  boundary  characteristic  function 


Now  putting 

Q(z,k)  ■  |  P(z.k)  1  + 1  R(z,k)  I, 

we  proved: 

Theorem  1.2  [20].  The  basic  approximation  (1 .8)  is  stable  if  Q{z,k)  >  0  for 
all 

{|z|  =  |  k |  »  1,  (z,k)  *  (1,1)}  u  0 < | k l <  1  >  .  (1.9) 

In  fact,  we  often  found  it  convenient  to  divide  the  (z.k)  domain  in  (1.9)  into 
three  disjoint  parts,  and  restate  Theorem  1 .2  as  follows: 

Theorem  1.2*.  Approximation  (1.8)  is  stable  if 


Q(z,k)>0  for  all  |z| 

=  j  K  |  =»  1,  K  *  1, 

(1.10a) 

Q(z,k»1)>0  for  ail 

|z|  -  1,  z  *  1, 

(1.10b) 

Q(z,k)  >0  for  all  |  z  | 

\t  1,  0  <  |  K  |  <  1 . 

(1.10c) 

The  advantage  of  this  setting  over  that  of  Theorem  1 .2  is  clarified  by  the 
following  lemma,  in  which  we  provide  helpful  sufficient  conditions  for  each  of 
the  three  inequalities  in  (1.10)  to  hold: 
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Lemma  1.3  [24]. 

(i)  Inequality  (1.10a)  holds  if  either  the  basic  scheme  (1.8a)  or  the 
boundary  conditions  (1 .8b)  are  dissipative. 

(ii)  Inequality  (1 . 1 0b)  holds  if  any  of  the  following  is  satisfied: 

(a)  The  basic  scheme  is  two-level. 

(b)  The  basic  scheme  is  three-level  and 

Q(z  =  -1,k  -  1)  >0.  (1.11) 

(c)  The  boundary  conditions  are  two-level  and  at  least  zero-order 
accurate  as  an  approximation  to  equation  (1.7). 

(d)  The  boundary  conditions  are  three-level,  at  least  zero-order 
accurate,  and  (1.11)  /s  satisfied. 

(iii)  Inequality  (1.10c)  holds  if  the  boundary  conditions  fulfill  the  von 
Neumann  condition,  and  are  either  explicit  or  satisfy 

Tide)  iJ^^O  V  0<|k|*1. 

i-o 


We  note  that  if  both  the  basic  scheme  and  the  boundary  conditions  are 
unitary  (i.e.,  strictly  nondissipative),  then  Q(z  ■  -1,k  »  -1)  ■  0;  hence  Theorem 
1 .2  is  rendered  useless.  For  such  cases  we  proved 

Theorem  1.4  [24].  Approximation  (1 .8)/'s  stable  if 

9P(z,k)  9P(z,k) 

dz  9k 


<0, 

Z  «  K  -  -1 


and  Q(z,k)  >  0  for  all 

{|z|-|k|-1,  (z,k)  #  ±(1,1)}U«(|Z|£1P  0<|K|<1} 

The  above  lemma  applies  to  this  theorem  precisely  in  the  same  way  it 
applied  to  Theorem  1 .2. 
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Tha  stability  criteria  obtained  in  Theorems  1.2  and  1.4  depend  both  on  the 
basic  difference  scheme  and  on  the  boundary  conditions,  but  not  on  the 
intricate  and  often  complicated  interaction  between  the  two.  Consequently, 
Theorems  1.2  and  1.4,  aided  by  Lemma  1.3,  provide  in  many  cases  a 
convenient  alternative  to  the  celebrated  stability  criteria  of  Kreiss  [31]  and  of 
Gustafsson,  Kreiss  and  SundstrCm  [26]. 

Having  the  new  criteria,  we  easily  established  stability  for  a  host  of  examples 
that  incorporate  and  generalize  most  of  the  cases  studied  in  recent  literature; 
e.g.,  [4,  5, 19,  20,  22-27,  30,  32,  38,  39,  42-  44,  47,  50].  To  mention  just  a  few  of 
our  examples,  we  proved  stability  for: 

(a)  Arbitrary  two-level  schemes,  with  boundary  conditions  generated  by 
either  the  explicit  or  implicit  one-sided  Euler  schemes. 

(b)  Arbitrary  two-level  schemes,  with  boundary  conditions  generated  by 
either  horizontal  extrapolation  or  by  the  one-sided  three-level  Euler 
scheme. 

(c)  Arbitrary  dissipative  schemes,  with  boundary  condition  generated  by 
oblique  extrapolation  or  by  the  Box  scheme. 

(d)  The  Crank-Nicolson,  Backward-Eufer,  Leap-Frog  and  Lax-Friedrichs 
schemes  (all  nondissipative),  with  boundary  conditions  generated  by 
either  oblique  extrapolative  or  by  the  one-sided  Weighted  Euler 
scheme. 


We  drew  great  satisfaction  from  the  fact  that  our  theory  and  examples  in  [19, 
20,  22-24]  were  used  already  by  a  number  of  authors,  including  Berger  [2], 
LeVeque  [34],  South,  Hafez  and  Gottlieb  [45],  Thun6  [49],  Trefethen  [50,  51], 
and  Yee  [53].  Thun6  [49],  in  his  effort  to  provide  a  software  package  for  stability 
analysis  of  finite  difference  approximations  to  hyperbolic  initial-boundary  value 
problems,  says:  "...Another  approach  has  been  to  derive  new  criteria,  based  on 
the  Gustafsson-Kreiss-SundstrSm  theory  but  more  convenient  for  practical  use... 
The  most  far-reaching  work  along  these  lines  has  been  made  by  Goldberg  and 
Tadmor  [19,  20, 221..." 

We  were  also  pleased  to  learn  that  part  of  our  theory  in  [24]  was  taught 
already  in  several  institutions  including  UCLA  and  the  University  of  Paris  VI. 


it 


2.  Multiplicativity  and  Stability  of  Matrix  and  Operator  Norms 

Let  V  be  a  normed  vector  space  over  the  complex  field  C,  and  let  2(V)  be  the 
algebra  of  bounded  linear  operators  on  V.  As  usual,  a  real-valued  function 

N  :  <8(V)  -*  R 

is  called  a  norm  on  s(V)  if  for  all  A,  B  e  ®(V)  and  a  e  C, 

N(A)  >0,  A*0, 

N{aA)  =  |  a  |  •  N(A), 

N(A  +  B)  £  N(A)  +  N(B). 

If  in  addition  N  is  multiplicative,  i.e., 


N(AB)  S  N(A)  N(B)  V  A,  B  6  5 B(V) , 

we  say  that  N  is  an  operator  norm  on  2(V).  If  ®(V)  is  an  algebra  of  (finite) 
matrices  and  N  is  multiplicative,  then  N  is  called  a  matrix  norm. 

The  first  multiplicative  example  that  comes  to  mind  is  of  course,  the  ordinary 
operator  norm 


j|  A  ||  «  sup  {  |  Ax  |  :  xe  V,  |  x  |  *  1  y. 


(2.1) 


where  |  •  |  is  the  vector  norm  on  V. 

If  V  is  a  (finite-  or  infinite-dimensional)  Hilbert  space,  then  perhaps  the  best 
known  example  of  a  nonmultiplicative  norm  on  <8(V)  is  the  numerical  radius 
(e.g.,  [1,6,  21, 28,  41]) 


r(A)  ■  sup  {  |  (Ax,  x)  |  :  x  e  V,  |  x  |  *  (x,  x) 


1/2 


-i> 


(2.2) 


12 


which  plays  an  important  role  in  stability  analysis  of  finite  difference  schemes  for 
multi-space-dimensional  hyperbolic  initial-value  problems  [21 , 33,  35,  52]. 


Another  example  of  considerable  interest  is  the  t  norm,  1  £  p  £  «»,  of  an 
n  x  n  complex  matrix  A  ■  (a  )  €  C 

r  j j  n  x  n 


(2.3) 


Ostrowski  [40]  has  shown  that  this  norm  is  multiplicative  (i.e.,  a  matrix  norm)  if 
and  only  if  1  £  p  £  2. 

Given  a  norm  N  on  «(V)  and  a  fixed  constant  p  >  0,  then  obviously  e  pN  is 
a  norm  too.  Clearly,  may  or  may  not  be  multiplicative.  -  If  it  is,  then  we  call  p  a 
multiplicativity  factor  lor  N.  That  is,  p  is  a  multi plicativity  factor  for  N  if  and  only  if 

N(AB)  £  pN(A)N(B)  V  A,  B  e  ®(V). 

Having  this  definition  one  can  obtain  at  once: 

Theorem  2.1  [36, 1 5].  Let  N  be  a  norm  on  s(V).  Then 
(i)  N  has  multiplicativity  factors  if  and  only  if 

p  i  *  Sup{N(AB)  :  N(A)  «  N(B)  ■  1 ;  A,  B  e  «(V) }  <  «».  (2.4) 


(ii)  If  p  .  <  oo,  then  p  is  a  multiplicativity  factor  for  N  if  and  only  /fp  £  p  .  . 

min  min 

In  the  finite-dimensional  case,  compactness  immediately  implies  that 
p  <  oo;  hence  N  always  has  multiplicativity  factors.  In  the  infinite- 

min 

dimensional  case,  however,  N  may  fail  to  have  multiplicativity  factors,  as  was 
demonstrated  by  Straus  and  myself  in  [1 51. 
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While  Theorem  2.1  seems  to  settle  the  question  of  characterizing 
multiplicativity  factors,  the  quantity  p  .  in  (2.4)  is  often  difficult  to  compute.  A 

more  practical  approach  towards  verifying  whether  a  constant  p  .  >  0  is  the 

min 

best  (least)  multiplicativity  factor  for  a  given  norm  N  is  implied  by  the  following 
obvious  observation: 

A  constant  p  .  >  0  is  the  best  (least)  multiplicativity  factor  for  N  if 


N(AB)  <;  p.  N(A)N(B)  V  A,  B  e  S(V), 


with  equality  for  some  nonzero  A  *  AQ,  B  =  BQ. 

With  this  observation  in  mind,  it  was  shown  by  Holbrook  [29]  (and 

independently  by  Straus  and  myself  in  [13])  that  if  V  is  a  Hilbert  space  of 

dimension  at  least  2,  and  if  r  is  the  numerical  radius  defined  in  (2.2),  then  pr  is 

an  operator  norm  on  «(V)  if  and  only  if  p  £  4;  i.e.,  the  best  multiplicativity  factor 

for  r  is  p  «  4. 
min 

Similarly,  Maitre  [36],  and  Straus  and  I  [17]  showed  that  the  best 
multiplicativity  factor  for  the  Cp  norm  on  Cn  x  n  defined  in  (2.3)  is 


mtn 


{nia* 


1  *  P  *  2, 
2£p£°°  . 


Often,  when  p  remains  unknown,  one  may  obtain  multiplicativity  factors 
min 

via  the  following  somewhat  stronger  version  of  a  result  by  Gastinel: 

Theorem  2.2  [3,  13].  Let  N  and  M  be  a  norm  and  an  operator  norm  on 
«(V),  respectively;  and  let  q  £  \  >  0  be  constants  such  that 

£M(A)  £  N(A)  £  qM(A)  V  A  e  S(V). 

2 

Then  any  p  with  p  £  r\!%  is  a  multiplicativity  factor  for  N. 
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This  result  was  utilized  by  Straus  and  myself  [13-16,  18]  to  obtain 
muitiplicativity  factors  for  certain  generalizations  of  the  numerical  radius,  called 
C-numerica!  radii. 

The  above  concepts  of  muitiplicativity  and  multiplicativity-factors  were 
extended  by  me  in  1983  as  follows: 

Definition  1.  Let  U,  V,  and  W  be  normed  vector  spaces  over  C;  and  let 
»  *  0(U,  W),  »  «(V,  W),  and  «3  ■  2(U,  V)  be  the  spaces  of  bounded  linear 

operators  from  U  into  W,  V  into  W,  and  U  into  V,  respectively.  If  Nv  N2,  and  N3 
are  norms  on  «1f  ®2,  and  ®3,  respectively,  and  p  >  0  is  a  constant  such  that 

N^AB)  <;  pN2(A)N3(B)  VAe  Be^, 

then  we  say  that  p  is  a  muitiplicativity  factor  for  N1  with  respect  to  N2  and  N3. 

In  analogy  with  Theorem  2.1  we  have  now: 

Theorem  2.3  [1 1].  Let  N2,  and  N3  be  norms  as  in  Definition  1 .  Then : 

(i)  N1  has  muitiplicativity  factors  with  respect  to  N2  and  N3  if  and  only  if 

pmjn  s  sup  {  Nt(AB)  :  N2(A)  =  N3(B)  *  1;  Ae  «2,  Be  $3}  <<*,. 

(ii)  If  [i  <  oo,  then  p  is  a  muitiplicativity  factor  for  N  with  respect  to  N 

min  *  ‘ 

and  N„  if  and  only  if  p  z  p  . 

3  min 

We  observe,  of  course,  that  a  constant  p  >  0  is  the  best  (least) 

min 

muitiplicativity  factor  for  N1  with  respect  to  N2  and  N3  if 

N^AB^p^AJN^B)  V  Ae  Sij.  Be  Sj, 
with  equality  for  some  nonzero  A  «  AQ,  B  -  BQ. 
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For  example,  if  V  is  a  Hilbert  space,  and  if  ||  •  ||  and  r  are  the  operator  norm 
and  numerical  radius  in  (2.1)  and  (2.2),  then  Holbrook  [29]  has  shown  that 


r(AB)  <;  2r(A)||  B  ||  V  A,  B  e  <8(V), 


with  equality  for  certain  A  ■  A.,  B  -  B„.  Thus,  u  «  2  is  the  best  multiplicativity 

O  0  min 

factor  for  r  with  respect  to  r  and  ||  •  j|. 


This  example  employs  only  a  single  vector  space  and  two  norms.  In  order  to 
demonstrate  the  idea  of  mixed  multiplicativity  to  its  full  extent,  consider,  for 
1  £  p  £  «»,  the  C  norm  of  an  m  x  n  matrix  A  ■  (a.)  e  C 

r  p  '  r  m  x  n 


Defining 


(2.5) 


I  proved: 

Theorem  2.4  [9].  Let  p,  q,  r  satisfy  1  £  p,  q,  r  £  <»,  and  let  q’  be  the 
conjugate  of  q  (i.e.,  1/q  +  1/q’  •  1).  Then  the  best  multiplicativity  factor  for  the 
C  norm  on  C  „  with  respect  to  the  C  norm  on  C  .  and  the  C  norm  on 

p  mxn  r  q  mxk  r 


mm 


Xp,(m)  y  n)  y<k). 
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That  is,  t or  all  A  e  C  .  and  B  e  C. 

mxk  kxn 

I  AB  Ip  £  XM(m)  y  n)  Xq.r(k)  |  A  |q  |  B  |r  ,  (2.6) 

where  this  inequality  is  sharp. 

Theorem  2.4  (which  generalizes  some  of  the  results  in  [7,8])  has  quite  a  few 
applications.  For  example  (see  [9, 12]),  taking  (2.6)  with  m  *  n  =  1,  we  get  an 
upper  bound  for  the  standard  inner  product  (x,  y)  on  Cn  in  terms  of  |  x  |p  and 
|  y  |q;  and  if  we  further  set  r  -  q‘  we  obtain  the  classical  Holder  inequality. 

Another  application  of  Theorem  2.4  concerns  the  "ordinary"  fp  operator-norm 


II A  llp 


xe  C 


(2.7) 


for  which  I  proved: 


Theorem 


A  € 


C  ..Be 

mxk 


2.5  [11]. 

Ckxn* 


Let  p,  q,  r  satisfy  1  £  p,  q,  r,  £  «». 


Then  for  all 


II  AB  ||P  <;  xM(m)  y(k)  y k)  y n)  n  a  nq  n  b  nr , 

where  the  inequality  is  sharp  if  either  q£p£r  or  r  £  p  £  q. 

Another  consequence  of  (2.6)  describes  the  equivalence  relation  between 
the  norms  in  (2.5)  and  (2.7): 

Theorem  2.6  [10].  Let  p,  q  satisfy  1  £  p,  q  £  and  let  q’  be  the  conjugate 
of  q.  Then  for  all  A  e  C 

^  mxn 
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I  A  Ip  ^pq(mn>l  A  Iq  * 

llAllp^Vm)Vn)  II  A  Hq  ’ 

llAllp^Vm)Wn),A|q- 

I A  |p  ^  (mn)1/ip|{  A  |(q  , 

where  the  first  three  inequalities  are  sharp. 

Having  the  above  results,  it  should  be  possible  to  construct  a  complete  table 
of  best  (least)  multiplicativity  factors  and  equivalence  constants  for  r,  |  •  |  and 
||  •  ||p,  as  well  as  for  other  useful  norms  such  as  the  Householder  norms 
described  in  [47]. 
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Memberships: 

Israel  Mathematical  Union 
American  Mathematical  Society 
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Grants  and  Awards: 

1 976-80  Principal  Investigator,  U.S.  Air  Force  Grant  AFOSR-76-3046 

1979-83  Principal  Investigator,  U.S.  Air  Force  Grant  AFOSR- 79-01 27 

1983-88  Principal  Investigator,  U.S.  Air  Force  Grant  AFOSR-83-0150 

1986-87  Distinguished  Lecturer  Award,  Technion 

1 988-  Principal  Investigator,  U.S.  Air  Force  Grant  AFOSR-88-01 75 

Publications:  See  attached  list. 

Talks  at  Universities  and  Research  Centers:  See  attached  list. 

Graduate  Students: 

Professor  Eitan  Tadmor,  M.Sc.,  1975 

Thesis:  "The  Numerical  Radius  and  Power  Boundedness" 

(Co-supervisor  with  Professor  G.  Zwas) 
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TALKS  AT  CONFERENCES  AND  MEETINGS 


Moshe  Goldberg 


1.  Invited  speaker,  international  Conference  on  Computational  Methods  in 
Nonlinear  Mechanics,  The  Texas  Institute  of  Computational  Mechanics, 
The  University  of  Texas  at  Austin,  Austin,  Texas,  September  1974,  title: 
"Stable  approximations  for  hyperbolic  systems  with  moving  boundary 
conditons". 

2.  Principal  speaker,  The  104th  Regular  Meeting  of  the  Association  for 
Computer  Machinery  (ACM),  Los  Angeles  Chapter,  Special  Interest  Group 
on  Numerical  Mathematics,  Los  Angeles,  California,  April  1975,  title: 
"Stable  approximations  for  hyperbolic  systems  with  moving  internal 
boundaries". 

3.  Invited  speaker,  American  Mathematical  Society  1975  Summer  Meeting, 
Special  Session  on  Numerical  Ranges,  Western  Michigan  University, 
Kalamazoo,  Michigan,  August  1975,  title:  "Inclusion  relations  between 
certain  sets  of  matrices". 

4.  Speaker,  The  1977  Army  Numerical  Analysis  and  Computer  Conference, 
Mathematics  Research  Center,  University  of  Wisconsin,  Madison, 
Wisconsin,  March  1977,  title:  "On  boundary  extrapolation  and  dissipative 
schemes  for  hyperbolic  problems". 

5.  Invited  speaker,  The  746th  American  Mathematical  Society  Meeting, 
Special  Session  on  Matrix  Theory,  California  State  University,  Hayward, 
California,  April  1977,  title:  "Some  inclusion  relations  for  c-numerical 
ranges". 

6.  Speaker,  The  1977  Dundee  Biennial  Conference  on  Numerical  Analysis, 
University  of  Dundee,  Dundee,  Scotland,  June  1977,  title:  "Dissipative 
schemes  for  hyperbolic  problems  and  boundary  extrapolation". 
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7.  Principal  speaker,  The  National  Science  Foundation  Conference  on 
Linear  and  Multilinear  Algebra,  University  of  California,  Santa  Barbara, 
California,  December  1977,  title:  "Numerical  ranges  and  numerical  radii”. 

8.  Speaker,  The  Eighth  U.S.  National  Congress  of  Applied  Mechanics, 
University  of  California,  Los  Angeles,  California,  June  1978,  title:  "Spectral 
analysis  of  hydroelastic  problems”  (with  R.S.  Chadwick). 

9.  Invited  speaker,  The  Second  International  Conference  on  General 
Inequalities,  Mathematical  Research  institute,  Oberwolfach,  West 
Germany,  August  1978,  title:  "Some  combinatorial  inequalities  and 
C-numerical  radii". 

10.  Principal  speaker,  Workshop  Series,  Five  One-Hour  Talks,  Institute  for  the 
Interdisciplinary  Applications  of  Algebra  and  Combinatorics,  University  of 
California,  Santa  Barbara,  California,  September  1979,  title:  "Numerical 
ranges  and  numerical  radii". 

11.  Principal  speaker,  The  October  1979  Meeting  of  the  Association  for 
Computing  Machinery  (ACM),  Los  Angeles  Chapter,  Special  Interest 
Group  in  Numerical  Mathematics,  Los  Angeles,  California,  October  1979, 
title:  "Stability  theory  for  difference  approximations  of  hyperbolic  partial 
differential  equations". 

12.  Invited  speaker,  The  1980  Annual  Meeting  of  the  Israeli  Society  for  the 
Applications  of  Mathematics,  Safad,  Israel,  May  1980,  title:  "Boundary- 
dependent  stability  criteria  for  difference  approximations  of  hyperbolic 
initial-boundary  value  problems". 

13.  Speaker,  The  1981  International  Conference  on  Convexity  and  Graph 
Theory,  University  of  Haifa,  Haifa,  Israel,  March  1981,  title:  "On  the 
convexity  of  numerical  ranges”. 

14.  Invited  speaker,  The  1981  Annual  Meeting  of  the  Israeli  Society  for 
Applications  of  Mathematics,  The  Weizmann  Institute  of  Science,  Rehovot, 
Israel,  April  1981,  title:  "Scheme-independent  stability  criteria  for 
difference  approximations  of  hyperbolic  initial-boundary  value  problems". 
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15.  Invited  speaker,  The  Third  International  Conference  on  General 
Inequalities,  Mathematics  Research  Institute,  Oberwolfach,  West  Germany, 
May  1981,  title:  "Better  stability  bounds  for  Lax-Wendroff  schemes  in 
several  space  dimensions". 

16.  Invited  speaker,  The  Toeplitz  Memorial  Conference,  Tel  Aviv  University, 
Tel  Aviv,  Israel,  May  1982,  title:  "The  numerical  radius:  from  Toeplitz  to 
modem  numerical  analysis”  (with  G.  Zwas). 

17.  Invited  speaker  (two  talks),  The  Fourth  International  Conference  on 
General  Inequalities,  Mathematics  Research  Institute,  Oberwolfach,  West 
Germany,  May  1983,  titles  (two  talks):  "New  inequalities  for  Cp  norms  of 
matrices",  and  "In  memoriam  Edwin  Beckenbach". 

18.  Invited  speaker,  The  AMS-SIAM  Summer  Seminar  on  Large-scale 
Computations  in  Fluid  Mechanics,  Scripps  Institute  of  Oceanography, 
University  of  California,  La  Jolla,  California,  June-July  1983,  title:  "New 
stability  criteria  for  difference  approximations  of  hyperbolic  initial-boundary 
value  problems”. 

19.  Invited  speaker,  The  1984  Annual  Meeting  of  the  Israel  Mathematical 
Union,  Applied  Mathematics  Session,  Tel  Aviv  University,  Tel  Aviv,  Israel, 
April  1984,  title:  "Convenient  stability  criteria  for  difference  approximations 
of  hyperbolic  initial-boundary  value  problems". 

20.  Speaker,  The  Gatlinburg  IX  Conference  on  Numerical  Algebra,  Univesity 
of  Waterloo,  Waterloo,  Canada,  July  1984,  titles  (two  talks): 
"Generalizations  of  the  Perron-Frobenius  Theorem",  and  "Norms  and 
multiplicativity". 

21.  Invited  attendee,  The  U.S.-lsrael  Binational  Workshop  on  the  Impact  of 
Supercomputers  on  the  Next  Decade  of  Computational  Fluid  Dynamics, 
Jerusalem,  Israel,  December  1984,  Panel  Discussion. 

22.  Invited  speaker,  The  1984  Haifa  Conference  on  Matrix  Theory,  Technion  - 
Israel  Institute  of  Technology  and  the  University  of  Haifa,  Haifa,  Israel, 
December  1984,  title:  "Submultiplicativity  of  matrix  norms  and  operator 
norms". 
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23.  Invited  speaker,  Joint  French-lsraeli  Mathematical  Symposium  on  Linear 
and  Nonlinear  Partial  Differential  Equations,  Numerical  Analysis,  and 
Geometry  of  Banach  Spaces,  The  Israel  Academy  of  Sciences  and 
Humanities,  Jerusalem,  Israel.  March  1985,  title:  "Stability  criteria  for 
difference  approximations  of  hyperbolic  initial-boundary  value  problems". 

24.  Principal  speaker,  Mathematics  Research  Conference,  California  Institute 
of  Technology,  Pasadena,  California,  October  1985,  title: 
"Submultiplicativity  of  matrix  norms  and  operator  norms". 

25.  Principal  speaker,  Southern  California  Functional  Analysis  Seminar 
(SCFAS),  California  State  University,  Los  Angeles,  California,  October 
1985,  title:  "Submultiplicativity  of  matrix  norms  and  operator  norms". 

26.  Invited  speaker,  The  1985  Haifa  Conference  on  Matrix  Theory,  Technion  - 
Israel  Institute  of  Technology  and  the  University  of  Haifa,  Haifa,  Israel, 
December  1985,  title:  "Submultiplicativity  and  mixed  submultiplicativity  of 
matrix  norms  and  operator  norms”. 

27.  Principal  speaker.  The  187th  Meeting  of  the  Association  for  Computing 
Machinery  (ACM),  Los  Angeles  Chapter,  Special  Interest  Group  in 
Numerical  Analysis,  Los  Angeles,  California,  February  1986,  title: 
"Stability  criteria  for  finite  difference  approximations  of  hyperbolic  initial¬ 
boundary  value  problems". 

28.  Invited  speaker,  The  Fifth  International  Conference  on  General 
Inequalities,  Mathematics  Research  Institute,  Oberwolfach,  West  Germany, 
May  1986,  title:  "Muttiplicativity  and  mixed-multiplicativity  of  operator 
norms  and  matrix  norms". 

29.  Speaker,  SIAM  Conference  on  Linear  Algebra  in  Signals,  Systems  and 
Control,  Boston,  Massachusetts,  August  1986,  title:  "Mixed  multiplicativity 
for  (p  norms  of  matrices". 

30.  Invited  speaker,  The  Third  Haifa  Matrix  Theory  Conference,  Technion  - 
Israel  Institute  of  Technology,  Haifa,  Israel,  January  1987,  title: 
"Equivalence  constants  for  f  norms  of  matrices". 
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31 .  Invited  speaker  and  Session  Chairman,  First  International  Conference  on 
Industrial  and  Applied  Mathematics,  Minisymposium  on  Finite  Difference 
Approximations  to  Hyperbolic  Initial-boundary  Value  Problems,  Paris, 
June-July  1987,  title:  "Convenient  stability  criteria  for  difference 
approximations  of  hyperbolic  initial-boundary  value  problems". 

32.  Invited  speaker,  Meeting  on  Numerical  Problems  for  Initial  and  Initial¬ 
boundary  Value  Problems,  Mathematics  Research  Institute,  Oberwolfach, 
West  Germany,  August  1987,  title:  "Stability  criteria  for  finite  difference 
approximations  to  hyperbolic  initial-boundary  value  problems". 

33.  Invited  speaker,  The  Fourth  Haifa  Matrix  Theory  Conference,  Technion  - 
Israel  Institute  of  Technology,  Haifa,  Israel,  January  1988,  title:  "On 
monotone  and  semi-monotone  matrix  functions". 

34.  Speaker  and  Session  Chairman,  Second  International  Conference  on 
Hyperbolic  Problems,  RWTH  Aachen,  Aachen,  West  Germany,  March 
1988,  title:  "Convenient  stability  criteria  for  difference  approximations  of 
hyperbolic  initial-boundary  value  problems". 

35.  Invited  speaker  and  Session  Chairman,  The  1988  Annual  Meeting  of  the 
Israel  Mathematical  Union,  Tel  Aviv  University,  Tel  Aviv,  Israel,  March 
1988,  title:  "Simple  stability  criteria  for  difference  approximations  of 
hyperbolic  initial-boundary  value  problems". 
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1.  T.J.  Watson  Research  Center,  IBM,  Yorktown  Heights,  New  York, 
Mathematics  Seminar,  September  1970. 

2.  NASA  Ames  Research  Center,  Moffet  Field,  California,  Thermo  and  Gas 
Dynamics  Divison,  Computation  Seminar,  November  1974. 

3.  University  of  California,  Berkeley,  California,  Department  of  Mathematics, 
Numerical  Analysis  and  Applied  Mathematics  Colloquium,  May  1975. 

4.  NASA  Langley  Research  Center,  Hampton,  Virginia,  Institute  for  Computer 
Applications  in  Science  and  Engineering  (ICASE),  ICASE  Seminar, 
August  1975. 

5.  University  of  California,  Santa  Barbara,  California,  Institute  for  the 
Interdisciplinary  Applications  of  Algebra  and  Combinatorics,  June  1976. 

6.  Tel  Aviv  University,  Tel  Aviv,  Israel,  Department  of  Mathematical  Sciences, 
Colloquium,  December  1976. 

7.  Tel  Aviv  University,  Tel  Aviv,  Israel,  Department  of  Mathematical  Sciences, 
Numerical  Analysis  Seminar,  December  1976. 

8.  The  Weizmann  Institute  of  Science,  Rehovot,  Israel,  Department  of 
Mathematics,  Colloquium,  December  1 976. 

9.  The  Hebrew  University,  Jerusalem,  Israel,  Institute  of  Mathematics, 
Colloquium,  December  1 976. 

10.  University  of  California,  Santa  Barbara,  California,  Department  of 
Mathematics,  Colloquium,  January  1 977. 

11.  University  of  California,  Los  Angeles,  California,  Department  of 
Mathematics,  Applied  Mathematics  Colloquium,  February  1977. 
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12.  Stanford  University,  Stanford,  California,  Computer  Science  Department, 
Numerical  Analysis  Seminar,  March  1977. 

13.  Case  Western  Reserve  University,  Cleveland,  Ohio,  Department  of 
Mathematics  and  Statistics,  Colloquium,  March  1977. 

14.  University  of  Southern  California,  Los  Angeles,  California,  Department  of 
Mathematics,  Colloquium,  November  1977. 

15.  University  of  California,  Los  Angeles,  California,  Department  of 
Mathematics,  Applied  Mathematics  Colloquium,  March  1978. 

16.  Polytechnic  Institute  of  New  York,  Brooklyn,  New  York,  Department  of 
Mathematics,  Colloquium,  March  1978. 

17.  Georgia  Institute  of  Technology,  Atlanta,  Georgia,  Department  of 
Mathematics,  Colloquium,  May  1978. 

18.  University  of  Georgia,  Athens,  Georgia,  Department  of  Mathematics, 
Colloquium,  May  1978. 

19.  Technion  -  Israel  Institute  of  Technology,  Haifa,  Israel,  Department  of 
Mathematics,  Colloquium,  October  1978. 

20.  Ben  Gurion  University,  Beer-Sheva,  Israel,  Department  of  Mathematics, 
Colloquium,  November  1978. 

21.  Tel  Aviv  University,  Tel  Aviv,  Israel,  Department  of  Mathematics, 
Colloquium,  November  1 978. 

22.  University  of  California,  Santa  Barbara,  California,  Institute  for  the 
Interdisciplinary  Applications  of  Algebra  and  Combinatorics,  Colloquium, 
May  1979. 

23.  Technion  -  Israel  Institute  of  Technology,  Haifa,  Israel,  Department  of 
Mathematics,  Analysis  Seminar,  three  one-hour  talks,  May  1980. 
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24.  University  of  California,  Santa  Barbara,  California,  Department  of 
Mathematics,  Linear  and  Multilinear  Algebra  Seminar,  six  75-minute  talks, 
August  1980. 

25.  Tel  Aviv  University,  Tel  Aviv,  Israel,  School  of  Mathematical  Sciences, 
Numerical  Analysis  Seminar,  May  1981. 

26.  Technion  -  Israel  Institute  of  Technology,  Haifa,  Israel,  Department  of 
Mathematics,  Analysis  Seminar,  June  1981. 

27.  University  of  California,  Santa  Barbara,  California,  Institute  for  the 
Interdisciplinary  Applications  of  Algebra  and  Combinatorics,  Colloquium, 
August  1981. 

28.  New  York  University,  New  York,  New  York,  Courant  Institute  of 
Mathematical  Sciences,  Numerical  Analysis  Seminar,  September  1981. 

29.  The  Hebrew  University,  Jerusalem,  Israel,  Institute  of  Mathematics, 
Colloquium,  February  1982. 

30.  The  Hebrew  University,  Jerusalem,  Israel,  Institute  of  Mathematics,  Partial 
Differential  Equations  Seminar,  March  1982. 

31 .  Tel  Aviv  University,  Tel  Aviv,  Israel,  Department  of  Mathematical  Sciences, 
Operator  Theory  Seminar,  April  1982. 

32.  University  of  California,  Santa  Barbara,  California,  Institute  for  the 
Interdisciplinary  Applications  of  Algebra  and  Combinatorics,  Colloquium, 
September  1982. 

33.  New  York  University,  New  York,  New  York,  Courant  Institute  of 
Mathematical  Sciences,  Numerical  Analysis  Seminar,  September  1983. 

34.  Tel  Aviv  University,  Tel  Aviv,  Israel,  Department  of  Mathematical  Sciences, 
Numerical  Analysis  Seminar,  January  1984. 

35.  Polytechnic  Institute  of  New  York,  Brooklyn,  New  York,  Department  of 
Mathematics,  Colloquium,  September  1984. 
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36.  University  of  California,  Los  Angeles,  California,  Department  of 
Mathematics,  Applied  Mathematics  Colloquium,  October  1985. 

37.  California  Institute  of  Technology,  Pasadena,  California, 

Department  of  Applied  Mathematics,  Colloquium,  November  1985. 

38.  California  Institute  of  Technology,  Pasadena,  California,  Department  of 
Mathematics,  Combinatorics  Seminar,  November  1985. 

39.  University  of  Southern  California,  Los  Angeles,  California,  Department  of 
Mathematics,  Colloquium,  March  1 986. 

40.  Centre  National  de  la  Recherche  Scientifique  et  University  Pierre  et  Marie 
Curie  (Paris  VI),  Paris,  France,  Numerical  Analysis  Seminar,  April  1986. 

41 .  £cole  Polytechnique,  Centre  de  Math6matiques  Apliqu6es,  Paris,  France, 
Applied  Mathematics  Seminar,  April  1986. 

42.  E=cole  Normale  Sup6rieure,  Paris,  France,  Centre  de  Math6matiques, 
Applied  Mathematics  Seminar,  April  1986. 

43.  Institute  National  de  Recherche  en  Informatique  et  en  Automatique  (INRIA), 
Rocquencourt,  France,  Seminar,  May  1 986. 

44.  University  of  Paris  IX,  Paris,  France,  Centre  de  Recherche  de 
Math6matiques  de  la  Decision,  Colloquium,  May  1 986. 

45.  New  York  University,  New  York,  New  York,  Courant  Institute  of 
Mathematical  Sciences,  Numerical  Analysis  Seminar,  September  1986. 

46.  Tel  Aviv  University,  Tel  Aviv,  Israel,  School  of  Mathematical  Sciences, 
Colloquium,  January  1987. 
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Theses: 

1 .  M.Sc.  Thesis,  "Quasi-conservative  hyperbolic  systems",  Tel  Aviv  University, 
Tel  Aviv,  Israel,  1970. 

2.  Ph.D.  Thesis,  "Stable  approximations  for  hyperbolic  systems  with  moving 
internal  boundary  conditions",  Tel  Aviv  University,  Tel  Aviv,  Israel,  1973. 

Published  Papers: 

1.  Numerical  solution  of  quasi-conservative  hyperbolic  systems  -  The 
cylindrical  shock  problem  (with  S.  Abarbanel),  Journal  of  Computational 
Physics  10  (1972),  1-21. 

2.  A  note  comparing  the  root  condition  and  the  resolvent  condition  (with  W.L. 
Miranker),  Information  Sciences  4  (1972),  285-288. 

3.  A  note  on  the  stability  of  an  iterative  finite-difference  method  for  hyperbolic 
systems,  Mathematics  of  Computation  27  (1973),  41-44. 

4.  Stable  approximations  for  hyperbolic  systems  with  moving  internal 
boundary  conditions  (with  S.  Abarbanel),  Mathematics  of  Computation  28 
(1974),  413-447. 

5.  Stable  schemes  for  hyperbolic  systems  with  moving  internal  boundaries 
(with  S.  Abarbanel),  in  "Computational  Mechanics",  edited  by  J.T.  Oden, 
Texas  Institute  for  Computational  Mechanics  (TICOM),  1974,  469-478. 

6.  On  matrices  having  equal  spectral  radius  and  spectral  norm  (with 
G.  Zwas),  Linear  Algebra  and  Its  Applications  8  (1974),  427-434. 

7.  The  numerical  radius  and  spectral  matrices  (with  E.  Tadmor  and  G.  Zwas), 
Linear  and  Multilinear  Algebra  2  (1975),  317-326. 
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8.  Numerical  radius  of  positive  matrices  (with  E.  Tadmor  and  G.  Zwas),  Linear 
Algebra  and  its  Applications  12  (1975),  209-214. 

9.  On  inscribed  circumscribed  conics  (with  G.  Zwas),  Elements  der 
Mathematik  31  (1976),  36-38. 

10.  Inclusion  relations  between  certain  sets  of  matrices  (with  G.  Zwas),  Linear 
and  Multilinear  Algebra  4  (1976),  55-60. 

11.  A  test  problem  for  numerical  schemes  for  nonlinear  hyperbolic  equations 
(with  S.  Abarbanel),  Computer  Methods  in  Applied  Mechanics  and 
Engineering  8  (1976),  331-334. 

12.  Inclusion  relations  involving  k-numerical  ranges  (with  E.G.  Straus),  Linear 
Algebra  and  Its  Applications  15  (1976),  261-270. 

13.  On  characterizations  and  integrals  of  generalized  numerical  ranges  (with 
E.G.  Straus),  Pacific  Journal  of  Mathematics  69  (1977),  45-54. 

14.  On  a  boundary  extrapolation  theorem  by  Kreiss,  Mathematics  of 
Computation  31  (1977),  469-477. 

15.  Elementary  inclusion  relations  for  generalized  numerical  ranges  (with  E.G. 
Straus),  Linear  Algebra  and  Its  Applications  18  (1977),  1-24. 

16.  On  a  theorem  by  Mirman  (with  E.G.  Straus),  Linear  and  Multilinear  Algebra 
5  (1977),  77-78. 

17.  On  boundary  extrapolation  and  dissipative  schemes  for  hyperbolic 
problems,  Proceedings  of  the  1977  U.S.  Army  Numerical  Analysis  and 
Computer  Conference,  ARO  Report  77-3  (1977),  157-164. 

18.  Scheme-independent  stability  criteria  for  difference  approximations  of 
hyperbolic  initial-boundary  value  problems.  1,  (with  E.  Tadmor), 
Mathematics  of  Computation  32  (1978),  1097-1107. 

19.  Norm  properties  of  C-numerical  radii  (with  E.G.  Straus)  Linear  Algebra 
and  Its  Applications  24  (1979),  113-132. 
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20.  On  certain  finite  dimensional  numerical  ranges  and  numerical  radii,  Linear 
and  Multilinear  Algebra  7  (1979),  329-342. 

21 .  Combinatorial  inequalities,  matrix  norms,  and  generalized  numerical  radii 
(with  E.G.  Straus),  in  "General  Inequalities  2",  edited  by  E.F.  Beckenbach, 
Birkhfiuser  Veriag,  Basel,  1980,  37-46. 

22.  Scheme-independent  stability  criteria  for  difference  approximations  of 
hyperbolic  initial-boundary  value  problems.  II,  (with  E.  Tadmor), 
Mathematics  of  Computation  36  (1981),  603-626. 

23.  On  the  numerical  radius  and  its  applications  (with  E.  Tadmor),  Linear 
Algebra  and  Its  Applications  42  (1 982),  263-284. 

24.  Operator  norms,  multiplicativity  factors,  and  C-numerical  radii  (with  E.G. 
Straus),  Linear  Algebra  and  its  Applications  43  (1982),  137-159. 

25.  On  the  mapping  A  ->  A+,  Linear  and  Multilinear  Algebra  12  (1983),285- 
289. 

26.  Combinatorial  inequalities,  matrix  norms,  and  generalized  numerical 
radii.  II,  (with  E.G.  Straus),  in  "General  Inequalities  3",  edited  by  E.F. 
Beckenbach  and  W.  Walter,  Birkhfiuser  Veriag,  Basel,  1983, 195-204. 

27.  Multiplicativity  of  £  norms  for  matrices  (with  E.G.  Straus),  Linear  Algebra 
and  its  Applications  52  (1983),  351-360. 

28.  Multiplicativity  factors  for  C-numerical  radii  (with  E.G.  Straus),  Linear 
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Research  of  Marvin  Marcus 
1983-1988 


L  Introduction 

This  report  covers  the  research  of  Marvin  Marcus  for  the  period  May  1 ,  1 983  - 
April  30,  1988  sponsored  by  the  Air  Force  Office  of  Scientific  Research,  grant 
number  AFOSR-83-0150. 

The  sequel  is  separated  into  the  following  sections: 

II.  General  Area  of  Research 

This  section  contains:  an  exposition  of  the  basic  mathematical  theory  of  the  finite  dimensional 
numerical  range;  a  description  of  two  algorithms  that  permit  effective  visualization  of  the 
structure  of  the  numerical  range;  an  example  of  the  implementation  of  algorithms  for 
visualizing  the  numerical  range  and  how  these  can  be  used  to  refute  or  substantiate  important 
conjectures;  a  list  of  continuing  problems  currently  under  investigation. 

III.  Research  of  M.  Marcus,  1983  - 1988 

This  section  contains  a  list  of  the  publications  completed  by  M.  Marcus  in  the  period  1983- 
1988  with  short  summaries  of  their  contents.  At  the  end  of  the  section  is  the  result  of  a 
computer  search  of  the  Science  Citation  Index  which  contains  the  total  number  of  references 
to  the  work  of  M.  Marcus  since  1983.  Self  references  have  been  excluded  in  the  search 
criteria.  This  data  provides  some  Quantitative  information  of  the  extent  to  which  the  research 
of  M.  Marcus  has  been  used  bv  other  investigators  working  in  the  general  area  of  appJied_and 
numerical  linear  algebra. 

IV.  Numerical  Range  Bibliography 

This  is  a  preliminary  version  of  a  bfoliography  of  779  citations  covering  the  numerical  range.  It 
has  been  sorted  alphabetically  by  first  author.  All  references  in  this  report  refer  to  the 
Bibliography.  We  are  currently  in  the  process  of  identifying  the  Mathematical  Reviews 
numbers  and  preparing  brief  summaries  of  each  of  the  papers.  In  view  of  the  very  large 
number  of  citations,  this  latter  project  will  probably  take  several  months  to  complete,  and  will 
be  incorporated  as  part  of  the  report  on  the  current  grant,  AFOSR-88-0175. 

V.  Appendix 

The  appendix  contains  the  vita  and  publication  list  of  Marvin  Marcus. 


B.  General  Area  of  Research 

Let  V  be  an  n -dimensional  unitary  space  and  let  A  be  a  (inear  transformation, 
A  :  V  V.  The  numerical  range,  or  field  of  values,  of  A  is  the  set  of  complex 
numbers 


W(A)  -  {(Ax.x)  |  l|  x  11  =  1>.  (1) 

The  numerical  radius  of  A,  w(A),  is  the  maximum  distance  of  any  point  in  W(A) 
from  the  origin.  By  choosing  an  orthonormal  (o.n.)  basis  of  V,  and  replacing  the 
inner  product  in  V  by  the  standard  inner  product  in  the  space  of  complex  column 
n-tuples,  the  computation  of  W(A)  is  reduced  to  an  equivalent  matrix  problem. 
Thus  we  assume  that  A  is  an  n-square  complex  matrix  and  that  the  inner 
product  of  two  column  n-tuples  (n-vectors)  is 

n 

(x.  y)  “  X  Vk*  (2) 

k- 1 


Elementary  results  concerning  W(A)  were  known  in  the  last  century  [165],  [398] 
and  in  the  first  decade  of  this  century  [69],  [317].  These  early  results  were 
usually  formulated  in  terms  of  bounding  rectangles  for  the  spectrum  of  A,  a( A), 
which  were,  in  fact,  containment  rectangles  for  W(A).  However,  it  was  not  until 
1918  and  1919  that  the  first  important  results  concerning  W(A)  were  proved  by 
Hausdorff  and  Toeplitz. 

it  is  a  classical  result  due  to  Hausdorff  [309]  and  Toeplitz  [712]  that  the 
numerical  range,  W(A),  is  a  convex  set.  Many  proofs  of  this  interesting  result 
have  appeared  in  the  intervening  years  since  the  original  Hausdorff-Toeplitz 
theorem  was  published.  Most  of  these  (e.g.,  see  [281])  depend  on  reducing  the 
problem  to  the  computation  of  the  numerical  range  of  a  2-square  matrix. 

There  have  been  a  number  of  interesting  papers  on  geometric  properties  of  the 
numerical  range  and  their  relation  to  the  similarity  invariants  of  A  (e.g.  [729], 
[41],  [175],  [213],  [351],  [376],  [728],  [561],  [144],  [154],  [184],  [49 7],  [583]).  From 
a  numerical  standpoint,  the  numerical  range  arises  in  many  contexts:  the 
constrained  eigenvalue  problem  [388];  the  theory  of  small  vibrations  [47],  [48]; 
Tchebychev  iteration  1  or  linear  systems  [446].  Much  of  the  interest  in  the 
numerical  range  of  a  matrix  A  is  motivated  by  the  fact  that  it  is  a  containment 
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region  for  the  spectrum  of  A.  In  fact,  for  normal  A,  W(A)  is  the  convex  hull  of  the 
spectrum  of  A.  It  might  be  conjectured  that  this  geometric  property  of  W(A)  is 
equivalent  to  A  being  normal.  In  fact,  M.  Marcus  and  B.N.  Moyls  [527]  showed 
that  for  n  s;  4  this  is  indeed  the  case,  but  for  n  >  4  it  is  not.  This  result  led  to  a 
sequence  of  related  papers  [19],  [216],  [29],  [30],  [463],  [474],  [486]  and  the 
introduction  of  a  class  of  operators  called  convexoid. 

The  numerical  range  of  any  linear  operator  is  the  union  of  the  numerical  ranges 
all  2-dimensional  real  compressions  of  A.  This  fact  is  the  basis  for  the  first 
algorithm  described  below.  If  1  £  k  £  n  and  P  is  a  k-dimensional  orthogonal 
projection,  then  the  restriction  of  PAP  to  the  range  of  P  is  called  a  k-dimensional 
compression  of  A.  For  k  *  2  and  A  an  n-square  complex  matrix,  a  2- 
dimensional  real  orthogonal  compression  of  A  is  the  2-square  matrix 

'  (Ax,  x)  (Av.x) " 

A  «  (3) 

™  L  (Ax- v)  (Av- v) . 

where  x  and  v  are  real  o.n.  column  n-tuples. 

The  following  are  well  known  properties  of  the  numerical  range.  The  set  W(A)  is 
unitarily  invariant  and  is  identical  with  the  set  of  all  diagonal  elements 
appearing  in  all  unitary  transforms  of  A  (i.e.,  in  all  matrices  unitarily  similar  to  A). 
The  numerical  range  of  every  principal  submatrix  of  A  is  a  subset  of  the 
numerical  range  of  A.  If  A  ■  B  ©  C  then  W(A)  *  H(W(B)  u  W(C)).  (H  denotes  the 
convex  hull.)  The  set  W(A)  is  a  closed  bounded  convex  region  of  the  plane 
containing  o(A),  the  spectrum  of  A,  i.e.,  containing  X  , ...,  X  ,  the  eigenvalues  of 

1  n 

•A.  Since  W(A)  is  convex,  it  also  contains 

P(A)  »  H(X1 . Xn).  (4) 

If  A  is  normal  then  W(A)  ■  P(A).  This  last  result  implies  that  if  A  is  normal  then 
the  extreme  points  of  W(A)  are  eigenvalues.  If  W(A)  «  {  X }  then  A  «  A.I  and  if 

W(A)  c  IRthen  A  *  A*,  i.e.,  A  is  hermitian. 

If  n  ■  2,  then  W(A)  is  an  ellipse  with  foci  the  eigenvalues  of  A;  if  A  has  the  form 


a 


*r 


0  Jtg 


then  the  length  of  the  semi-minor  axis  of  the  ellipse  is  |  a  |/2.  The  precise 
equation  for  the  boundary  of  the  numerical  range  of  a  non-normal  matrix  for 
n^2  has  been  given  by  Murnaghan  [528]  and,  in  more  explicit  form,  by 
Kippenhahn  [376].  Kippenhahn  also  gives  bounds  for  the  diameter  of  W(A).  M. 
Feidler  obtained  [213]  an  equation  for  the  boundary  of  W(A). 

Since  the  numerical  range  contains  the  spectrum  of  A  it  is  of  considerable 
importance  from  the  standpoint  of  eigenvalue  localization.  In  fact,  this  was  the 
starting  point  for  a  number  of  papers  on  classical  eigenvalue  localization  theory 
including  work  by  W.V.  Parker  [554],  [556]  and  A.B.  Famell  [206],  [208]. 

P.  Henrici  [312]  related  the  distance  between  the  boundary  of  W(A)  and  P(A) 
with  a  measure  of  the  departure  of  A  from  normality. 


In  a  paper  written  in  1952  [254]  W.  Givens  defined  for  As  Mn(<D)  the  following 
set: 


Fh(A) 


f  (HAx,  x) 
\  (Hx,  x) 


X6Cn|, 


H  p.d. 


Givens  proved  that  if  H  is  p.d.  (positive  define  hermitian),  H  *  T*T,  then  FU(A)  = 

W(TAT  "1).  He  also  showed  that  if  A  has  an  elementary  divisor  of  degree  at  least 
2  associated  with  the  root  X,  then  X  is  an  interior  point  of  FH(A)  for  every  p.d.  H. 

An  immediate  consequence  of  this  last  result  is  that  if  X  is  an  eigenvalue  on  the 
boundary  of  W(A)  then  X  occurs  only  in  linear  elementary  divisors.  Givens’ 
main  result  was 


P(A)  .  O  Fh(A). 

H  p.d. 

He  also  showed  that  a  necessary  and  sufficient  condition  that  FU(A)  =  P(A)  for 

n 
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some  p.d.  H  is  that  the  elementary  divisors  corresponding  to  roots  on  the 
boundary  of  P(A)  are  all  linear.  Givens'  results  suggest  that  for  an  appropriate 
choice  of  H  one  might  obtain  information  about  the  eigenvalues  of  A, 
particularly  those  on  the  boundary  of  P(A),  from  properties  of  the  numerical 
range. 

In  [175]  W.F.  Donoghue  proved  that  every  non-differentiable  boundary  point  of 
W(A)  is  an  eigenvalue  of  A. 

0.  Taussky  [690]  has  shown  that  if  A  *  0  and  tr(A)  »  0  then  0  is  in  the  interior  of 
W(A). 

A  result  of  C.R.  Putnam  [580]  states  that  if  C  »  AB  -  BA  then  0  is  in  the  interior  of 
W(C). 

Ballantine  [41]  has  presented  a  series  of  algorithms  to  determine  for  a  given 
complex  number  z  and  a  given  A  e  Mn(C)  whether  or  not  z  is  in  W(A),  z  is  a 
boundary  point  of  W(A),  or  z  is  an  extreme  point  of  W(A). 

In  a  paper  written  in  1963  [457]  M.  Marcus  and  R.C.  Thompson  examined  the 
numerical  range  of  the  Hadamard  product  A  *  B  of  two  matrices.  They  showed 
that  if  A  and  B  are  normal  and  c^,  ...,  and  j^, ....  (5^  are  the  eigenvalues  of  A 

and  B  respectively  then  W(A  *  B)  is  a  subset  of  the  convex  polygon  spanned  by 
(a  0.  +  a  | (J  )/2,  1  £  i  £  j  £  n.  This  result  was  used  to  yield  localization  theorems 
for  permanents  and  determinants. 

T.  Saitd  [600]  considered  the  question:  When  is  the  relation  W(A  ®  B)  * 
H(W{A)W(B))  valid?  In  a  particular  answer  to  this  question  he  proved  that  if 
W(A  ®  B)  ■  H(o(A  ®  B))  then  the  above  equality  holds.  He  also  showed  that  in 
general  H(W(A)W(B))  c,  W(A  ®  B)  but  that  there  exist  A  and  B  for  which  the 
inclusion  is  strict. 

in  a  series  of  papers  [348],  [349],  [351  A],  [356]  Johnson  examined  various 
inclusion  relations  involving  W(A).  In  [351],  for  n  ■  2,  he  determined  the  major 
and  minor  axes  of  the  ellipse  W(A)  in  terms  of  the  entries  of  A  when  A  is  real. 
He  then  utilized  this  result  to  determine 
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S(A)  -  H(UW(A[i,j|i,j]) 

*»J 


for  A  n-square  real  and  showed  that  S(A)  c  W(A). 

In  [354],  [353]  Johnson  studied  the  Hadamard  product  A  *  B  of  A  and  B(A*B« 

[a^).  He  proved  that  if  A  e  Mn(C)  and  for  some  0  £  0  £  2ic,  e^H  is  p.d.  then 

W(H  ®  A)  «  W(H)W(A).  Furthermore,  if  A  e  Mn(C)  and  N  e  Mm(C)  is  normal 

then  W(N  ®  A)  c  H(W(N)W(A)).  A  corollary:  if  N  and  A  are  in  Mn(£)  and  N  is 

normal,  then  W(N  *  A)  c  H(W(N)W(A));  if  further  N  is  p.d.  then  W(N  •  A)  c 
W(N)W(A). 

Since  the  effective  visualization  of  the  set  W(A)  has  been  an  important  part  of 
this  research  it  is  important  to  be  aware  that  several  algorithms  exist  for 
graphing  the  convex  hull  of  a  set  of  points.  Sedgewick  describes  the 
implementation  of  the  package  wrapping  algorithm  in  [Algorithms,  2nd  ed., 
Robert  Sedgewick,  Addison  Wesley,  1988]  which  is  not  unrelated  to  one  of  the 
algorithms  developed  below  to  visualize  the  numerical  range: 

1 .  Find  the  point  with  the  least  y  coordinate. 

2.  Imagine  a  horizontal  line  through  this  point. 

3.  Sweep  that  horizontal  line  through  a  positive  angle  Q  until  it  intersects  with 
another  point. 

4.  Add  that  point  to  the  boundary  of  the  convex  hull. 

5.  If  the  new  point  is  not  the  starting  point,  goto  step  2. 

► 

Obviously  this  algorithm  is  suitable  for  finite  sets  of  points  only. 

To  visualize  the  numerical  range  it  is  required  to  graph  its  boundary.  In  the 
second  algorithm  described  below,  an  effective  means  of  computating  the 
boundary  of  W(A)  is  described. 

The  mathematical  results  at  the  basis  for  the  visualization  algorithms  will  be 
discussed  next.  Let  A  e  Mn((D),  let  B  be  any  principal  submatrix  of  A  and  let  U 

€  Mn(<D)  be  any  unitary  matrix.  Also  let  c(A)  denote  the  spectrum  of  A,  i.e.,  a(A) 
is  the  set  of  all  eigenvalues  of  A.  Then 
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W(cA)  -  cW(A), 

(5) 

W(cln+A)  -  c  +  W(A), 

(6) 

W(A*)  -  W(A), 

(7) 

W(B)  c  W(A), 

(8) 

W(U*AU)  -  W(A), 

0) 

W(A)  cIR  iff  A  is  hermitian, 

(10) 

W(A)  c  ilR  iff  A  is  skew-hermitian, 

(11) 

W(A) « { 0 }  iff  A  =  0, 

(12) 

c(A)  c  W(A), 

(13) 

W(A)  »  { c} iff  A  =  cln. 

(14) 

The  following  theorem,  known  as  the  elliptical  range  theorem,  completely 
describes  the  structure  of  W{A)  for  A  e  M2(C).  It  is  the  basis  for  proving  the  fact 
that  W(A)  is  always  convex,  the  most  important  theorem  about  the  numerical 
range.  It  is  also  the  basis  of  an  effective  visualization  algorithm. 

Theorem  1.  Let  A  e  M2(C)  with  eigenvalues  X  and  p.  Define  the  following 
numbers  associated  with  A: 


1/2 


Ui  - 1 


I  aii  I 


,2  .  .  ,2  .  .2,1/2 

a  a  (v  -jX|  - 1  p  |  )  . 


(15) 

(16) 


Then  W(A)  is  an  elliptical  region  bounded  by  an  ellipse  (possibly  degenerate) 
whose  description  is  as  follows: 


foci:  Xp; 

.  .  .  (v2-2ReXp)1/2 

semi-major  axis:  - - —  — — 

semi-minor  axis:  ^ 


(17) 

(18) 
(19) 


TTC- 


v 
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An  important  result  proved  in  [489]  is  contained  in 


Theorem  2.  Let  A  be  an  n-square  complex  matrix.  Then  W(A)  is  the  union  of  all 
the  sets 


W(AJ 


(20) 


as  x  and  v  run  over  all  pairs  of  o.n.  vectors  in  Mn  1  (R). 

From  a  computational  standpoint  it  is  important  to  note  that  although  W(A) 
consists  of  complex  numbers  of  the  form  y*Ay,  ye  M_  1  (C),  it  is  nevertheless 

n*  1 

the  case  that  only  isal  x  and  v  are  required  in  (20). 

We  can  state  the  following  useful  result  for  computing  the  numerical  radius  of  a 
2x2  matrix. 


Theorem  3.  Let  A  be  unitarily  similar  to 


X  a 

0  ii 


a£0. 


(21) 


For  s  e  [0, 1]  define  the  function 

d(s)  ■  |  sX  +  (1  -  s)p  |  +  a  Vs(1  -  s)  (22) 


then 


w(A)  -  max  d(s)  (23) 

where  the  max  in  (23)  is  computed  for  s  e  [0, 1]. 


For  matrices  A  e  M2(C)  which  are  unitarily  similar  to  a  real  matrix  it  is  possible 
to  give  an  explicit  formula  for  w(A)  in  terms  of  the  entries  of  A.  In  the  following 
theorem,  w(A)  is  explicitly  exhibited  in  terms  of  the  entries  in  the  upper 


triangular  form  of  A  for  a  2-square  matrix.  This  result  is  useful  for  determining 
an  approximation  from  below  for  the  numerical  radius  of  an  arbitrary  A. 


Theorem  4.  If  A  e  M2(<C),  upper  triangular, 


A 


X  a 
0  H 


a£0. 


and  if  A  is  unitarily  similar  to  a  real  matrix  then  the  numerical  radius  w(A)  can  be 
determined  as  follows: 

I.  X  and  p  are  real.  Then  w(A)  is  the  larger  of  the  two  numbers 


|X  +  p±V (X -  p)2  +  a2 1  . 

2 

II.  X  and  p  are  complex  conjugates:  X  »  h  +  ik.  p  *  h  -  ik,  k  *  0.  If  2k2  £  a|  h  |  then 


w(A) 


JM 

|k| 


JJl 


4k 


III.  X  and  p  are  complex  conjugates:  X  »  h  +  ik,  p  *  h  -  ik,  k  *  0.  If  2k2  <  oc|  h  | 
then 


w(A)-|hj+£. 

If  C  is  a  convex  set  in  C  which  is  closed,  i.e.,  contains  all  its  limit  points,  and 
which  is  not  all  of  <D,  then  C  is  the  intersection  of  all  its  supporting  half-planes. 
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This  fact  is  geometrically  evident  and  it  is  not  difficult  to  prove.  The  geometry  of 
this  situation  is  very  useful  in  developing  an  algorithm  for  constructing  W(A). 
Theorem  2  provides  us  with  one  effective  method  for  constructing  W(A).  The 
present  discussion  will  lead  us  to  another  such  method.  Thus,  let  W(A)  be  the 
numerical  range  of  an  arbitrary  A  e  Mn(C). 


(24) 


The  idea  is  simple:  we  want  to  construct  a  relatively  dense  set  of  support  lines 
for  W(A).  Then  W(A)  will  be  accurately  depicted  as  the  intersection  of  the 
corresponding  support  half-planes.  In  fact,  simply  drawing  a  sufficiently  dense 
set  of  such  support  lines  will  define  W(A)  with  great  accuracy. 


r* - 

I 

I 


(25) 


Of  course,  the  problem  is  to  devise  a  computationally  reasonable  method  of 
determining  the  support  lines  L.  We  will  make  our  method  depend  on 
computing  the  dominant  eigenvalue  of  a  sequence  of  appropriate  hermitian 
matrices.  Thus  let  L  be  a  fixed  but  arbitrary  support  line  for  W(A).  Then  perform 
a  counterclockwise  rotation  in  the  plane  through  an  angle  e  chosen  so  that  the 
rotated  image  of  L,  call  it  Le,  is  perpendicular  to  the  x-axis  (see  (24)).  Clearly 

each  such  support  line  L  determines  a  unique  Lfl.  For  a  given  6,  if  we  can 
determine  the  equation  for  L.,  then  the  equation  for  L  is  obtained  by  elementary 

geometry.  Now,  Lfl  is  a  support  line  for  W(e,0A)  =  el0W(A).  Write  A(0)  =  e^A  and 
let  H(8)  and  K(9)  be  the  hermitian  parts  of  A(0): 

A(0)  »  H(0)  +  iK(9). 

Then  if  u*u  ■  1 , 

u*A(8)u  =  u*H(0)u  +  iu*K(0)u 
and 
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Re  u*A(0)u  =  u*H(0)u. 

Hence 

max  Reu*A(0)u  *  max  u*H(0)u. 

J*U  >1  u*u  -  1 

But  H(0)  is  hermitian  and  thus 

max  u*H(0)u  =  A,(0) 

u*u  -  1 

where  A.(0)  is  the  largest  eigenvalue  of  H(0).  In  fact,  the  maximizing  u  is  an 
eigenvector  of  H(0)  corresponding  to  X(0).  Assume  that  it  is  feasible  to  compute 
X(0).  Then  the  equation  of  L.  is  obviously 

x  =  X(0), 

or  in  complex  number  notation 

Rez  a  X(0). 

Now  a  point  z  =  |  z  le^  *  x  +  iy  lies  on  the  line  L  iff  e^z  lies  on  the  line  La,  i.e,  iff 

Ree^z  =  X(0), 

Re  e®!  z  |ei<p  ■  X(0), 

Re|z|e'(e+<p)  =  X(0), 

|z|cos(0  +  <p)  »  X,(0), 

|  z  |  costp  cos0  - 1  z  |  sin<p  sin0  *  X(0), 

xcos0-ysin0  -  X(0).  (26) 

Thus  (26)  is  the  equation  for  L  in  rectangular  coordinates.  The  line  (26)  is 
known  once  0  is  specified  and  \{B)  is  computed.  A  sensible  scheme  might  be 


53 


to  specify  a  sequence  of  N  values  of  6  of  the  form 


e*  =  k7P  k-0'1 . N-' 

and  construct  the  lines 

Lk  :  x  cosek-ysin9k  *  X(9k),  k  =  0 . N  - 1. 

The  method  for  depicting  W(A)  just  described  can  also  be  used  to  determine  the 
numerical  radius  w(A).  In  fact,  we  can  easily  verify  that 

w(A)  «  max  X(0). 

0  e  [0,  2n] 


Algorithm  1 

The  first  algorithm  is  based  on  Theorem  1,  Theorem  2,  and  Theorem  4. 
Theorem  1  is  the  elliptical  range  theorem.  Theorem  2  states  that  W(A)  is  the 
union  of  the  numerical  ranges  of  all  the  real  2-square  orthogonal  compressions 
of  A.  Theorem  4  is  the  explicit  formula  for  w(A)  of  a  2-square  matrix  unitarily 
similar  to  a  real  matrix.  The  algorithm  approximates  W(A)  from  the  inside. 

1  Generate  a  random  o.n.  pair  of  vectors,  x  and  v 

2.  Compute  Axy 

3.  Apply  the  elliptical  range  theorem  to  Aw  to  obtain  W(Axy) 

4.  Graph  WjA^) 

5.  Update  w(A)  with  maximal  value  of  wjA^) 

6.  Goto  1 . 

The  Macintosh  has  many  built  in  ROM  routines  for  drawing  objects  on  the 
screen.  The  routines  that  were  used  in  this  program  were  paintoval  and  lineto. 
The  paintoval  command  draws  an  oval  inside  a  specified  rectangle.  This 
rectangle  is  situated  in  the  plane  with  its  sides  parallel  to  the  axes.  Thus,  it  was 
not  possible  to  draw  an  inclined  ellipse.  The  foci  of  the  ellipse  had  to  lie  along 
the  real  or  the  imaginary  axis.  This  means  that  the  eigenvalues  of  the  2-square 
compressions  of  A  had  to  be  real  or  complex  conjugates  of  one  another.  This 
fact  restricted  our  ability  to  quickly  depict  W(A)  for  an  arbitrary  complex  matrix 
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using  this  algorithm. 


Step  1  consists  of  generating  a  pair  of  o.n.  vectors.  This  is  done  by  randomly 
generating  two  n-vectors  with  components  e  [-1,  1].  These  two  vectors  are 
checked  to  make  sure  their  lengths  are  greater  than  0  and  then  orthonormalized 
using  the  Gram-Schmidt  process. 

Step  2  consists  of  computing  A^.  This  is  done  by  computing 

(Ax,x)  (Av,x) ' 

(Av,x)  (Av.v) 

When  it  was  necessary  to  compute  a  vector-matrix-vector  product  as  above,  the 
operations  were  applied  as  follows: 

(Ax,x)  *  (x*(Ax)}. 

Step  3  consists  of  applying  the  elliptical  range  theorem  to  Axy.  This  entailed 
computing  the  eigenvalues  of  A^ ,  rl  and  r2.  The  eigenvalues  are  computed  by 
solving  the  characteristic  polynomial  of  Axy.  These  eigenvalues  are  the  foci  of 
the  ellipse  which  is  the  numerical  range  of  Axy.  Next,  the  value  alpha  is 
computed.  This  is  the  length  of  the  minor  axis.  Depending  on  the  values  of 
alpha,  rl,  and  r2,  W(Axv)  has  different  features  and  is  graphed  accordingly. 

Step  4  graphing  W(Axy).  If  alpha  ■  0  then  W(Axv)  collapses  to  a  line  segment 
joining  rl  and  r2.  If  this  is  the  case  is  drawn  using  a  straight  line. 

If  rl  and  r2  are  complex  conjugates,  alpha  >  0,  then  W(Axv)  is  situated  in  the 
plane  with  its  major  axis  parallel  to  the  imaginary  axis. 

If  rl  an  r2  are  real,  alpha  >  0,  then  W(Axv)  is  situated  in  the  plane  with  its  major 
axis  along  the  real  axis. 

After  each  W(Axy)  is  drawn,  it  is  checked  to  see  whether  any  of  its  points  are 
either  the  topmost,  bottommost,  leftmost,  or  rightmost  points  in  W(A)  exhibited  so 
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far.  After  each  iteration,  the  convex  polygon  is  drawn  connecting  these  four 
extreme  points.  This  speeds  the  approximation  of  the  convex  hull  of  A  from 
inside. 


To  illustrate  this,  suppose  we  are  trying  to  approximate  the  numerical  range  of 
the  matrix 


A 


0  1  0 
0  0  1 
1  0  0 


We  know  that  W(A)  is  a  triangle  (A  is  normal  and  its  numerical  range  is  the 
convex  hull  of  its  eigenvalues).  After  two  iterations  we  may  have  two  ellipses 
situated  as  depicted  below,  that  approximate  the  triangle  from  the  inside.  If  we 
connect  the  extreme  points  on  the  ellipses  we  get  a  closer  approximation  to 
W(A).  Joining  the  extreme  values  is  performed  at  every  iteration  past  the  first 
one. 


Step  5.  update  w(A).  Theorem  4  provides  us  with  a  closed  form  formula  for 
evaluating  w(A)  for  a  2-square  matrix  unitarily  similar  to  a  real  matrix.  The 
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theorem  has  three  alternatives:  I)  the  eigenvalues  are  real;  II)  the  eigenvalues 
are  complex  conjugates  of  one  another,  h  ±  ik  and  2k2  £  a  |  h  |;  and  111)  the 
eigenvalues  are  complex  conjugates,  h  ±  ik,  and  2k2  <  a  |  h  |.  The  conditions  of 
the  theorem  are  checked  against  the  eigenvalues,  rl  and  r2,  and  alpha  to  see 
which  case  holds.  Then  the  value  of  w(Axv)  is  computed.  This  is  compared  to 
the  maximum  value  to  date,  and  the  maximum  value  is  updated  if  necessary. 

Step  6  -  Goto  Step  1.  This  program  has  no  set  stopping  criteria.  It  is 
programmed  to  run  indefinitely.  Theorem  2  states  that  W(A)  is  the  union  of  all 
WfA^)  and  x  and  v  are  being  generated  randomly.  When  the  image  of  W(A) 

appears  to  have  stabilized  into  a  convex  shape  then  it  is  interrupted. 

Algorithm  2 

The  second  algorithm  is  based  the  fact  that  the  numerical  range  is  a  convex  set. 
It  implements  the  algorithm  that  visualizes  W(A)  as  the  intersection  of  half¬ 
spaces  of  support  lines.  The  algorithm  approximates  W(A)  from  the  outside. 

The  algorithm  goes  as  follows: 

1 .  Determine  an  angle  y 

2.  n  <-  trunc  (2n/y  +  0.5) 

3.  for  j :»  0  to  n  do 

1.  0«-j»y 

2.  + 

3.  w  <-  X  (H{6)) 

max 

4.  maxw  <-  max(w,  maxw) 

5.  graph  the  support  line  corresponding  to  X  (H(6)) 

max 

Unlike  the  implementation  of  Algorithm  1 ,  this  algorithm  is  able  to  exhibit  W(A) 
for  an  arbitrary  complex  matrix.  It  is  not  restricted  to  matrices  that  are  unitarily 
similar  to  a  real  matrix. 

Step  1 .  The  user  is  able  to  enter  any  choice  for  y  e  [0,  2 tcJ. 

Step  2.  Here  we  compute  the  number  of  iterations  for  the  program.  This  is 
unlike  the  first  algorithm.  This  program  will  terminate  after  a  predetermined 
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number  of  iterations  with  an  outer  approximation  to  W(A)  and  an  upper  bound 
for  w(A).  The  smaller  the  angle  y  specified,  the  more  iterations  and  the  better 
the  approximation  to  W(A). 

Step  3.1.  Each  successive  angle  0  is  computed. 

Step  3.2.  H(0),  the  hermitian  part  of  e®A,  is  calculated. 


Step  3.3.  The  power  method  is  run  on  H(0)  to  determine  its  largest  eigenvalue. 
The  version  of  the  power  method  implemented  here  is  the  Rayleigh  Quotient 
method.  This  method  will  find  the  largest  eigenvalue  in  modulus  of  the  given 
matrix.  For  this  algorithm  the  rightmost  eigenvalue  is  required.  To  get  around 
this  problem  the  matrix  T  was  computed,  T  =  H(0)  +  [|  H(0)  ||^  K  (Here  ||  A  ||  is 


the  1  -norm,  max  ^  |  aa  I  i  =  1 . n.)  This  ensures  that  T  is  a  positive  semi- 

‘j-i 

definite-matrix  (X.  2:  0,  i  »  1,  ....  n).  Thus  the  rightmost  eigenvalue  of  T  is  the 

eigenvalue  of  maximum  modulus.  The  rightmost  eigenvalue  of  H(0)  was 
computed  by  X  (H(0))  =  X  (T)  -  II  H(0)  ||  . 

max  ma*  1 


max 


Step  3.4.  The  maximum  of  the  values  X  (H(0)),  0  e  [0,  2jc],  is  an 

max 

approximation  to  the  numerical  radius  of  A.  This  value  is  maximized  at  every 
iteration. 


Step  3.5.  The  equation  of  the  support  line  at  the  point  e',eX  (H(0)),  rotated 

max 

counterclockwise  through  0,  is  x  »  Xmax(H(0)).  Thus  the  equation  of  the  support 
line  itself  is 


xcos0-ysin0  =  X(H(0)). 


A  problem  was  encountered  with  the  implementation  of  this  algorithm.  In  Step 
3.3  when  the  power  method  is  applied,  an  initial  estimation,  xQ,  of  an 

eigenvector  is  required.  Originally,  the  code  was  written  so  that  xQ  =  [1, ....  1]T. 

This  presented  no  problem  with  the  majority  of  examples  for  which  the  algorithm 
was  tested.  But  the  program  consistently  failed  for  any  doubly  stochastic  matrix. 
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This  failure  can  be  explained  however,  the  simplest  solution  computationally  is 
to  generate  a  random  starting  vector. 

Nilpotent  Matrices 

Using  the  visualization  algorithms  it  is  easy  to  construct  examples  of  nilpotent 
matrices  whose  numerical  ranges  are  disks  centered  at  the  origin.  The 
question  arises:  what  are  necessary  and  sufficient  conditions  on  a  nilpotent  A 
so  that  W(A)  is  a  disk  centered  at  the  origin? 

Let  A  to  be  an  arbitrary  n-square  matrix.  There  is  no  loss  in  generality  in 
assuming  that  0  e  W(A). 

Theorem  5.  Let  A  «  H  +  i  K  be  the  hermitian  decomposition  of  A.  Let  \(Q)  be 
the  maximum  eigenvalue  of 


cos  6  H  -  sin  0  K  (27) 

If  0  €  W(A)  then  W(A)  is  a  disk  centered  at  the  origin  iff  the  maximum 
eigenvalue  X(0)  of  (27)  is  independent  of  0,  0  £  0  £  2k. 

Theorem  6.  Let  A  be  an  n-square  real  nilpotent  matrix.  For  n  =  3,  W(A)  is  a  disk 
centered  at  the  origin  iff 


1r((A2)TA)  =  0. 

For  n  ■  4,  W(A)  is  a  disk  centered  at  the  origin  iff 

tr<(A*)TA)  -  0 


and 


tr((A3)TA)  =  0. 

Research  is  currently  underway  to  extend  Theorem  6  to  general  n-square 
matrices. 
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Normal  Matrices  and  Symmetry 


Let  A  be  a  linear  operator  on  a  finite  dimensional  unitary  space  V  of  dimension 
n.  The  kth  higher  numerical  range  of  A,  denoted  by  Wk(A),  is  the  totality  of 

complex  numbers  tr(PAP)  where  P  runs  over  all  k-dimensional  orthogonal 
projections  on  V.  Very  recently  [490]  we  were  able  to  prove  that  Wk(A)  is 
polygon  with  the  real  axis  as  a  line  of  symmetry,  k  «  1 ,  n,  if  and  only  if  A  is 
normal  with  a  real  characteristic  polynomial.  We  also  constructed  several 
nonnormal  examples  to  investigate  the  extent  to  which  the  symmetry  of  all  the 
WR(A)  is  required. 


Visualization  of  the  Numerical  Range 

As  an  example  of  the  use  of  the  visualization  algorithms  described  above, 
consider  the  following  problem.  Is  it  the  case  that 

Wk{A)«Wk(B),  k-1 . n  (28) 

suffice  to  conclude  that  the  two  n-square  matrices  A  and  B  are  unitarily  similar? 
Recall  that 


W.  (A) 


i 

j-i 


(Ax.,  xp, 


xko.n. 


so  that  W^A)  is  simply  W(A).  To  investigate  this  conjecture  we  take  n  *  3  and 
then  directly  confirm  that  the  conditions  (28)  are  equivalent  to 


W(A)  -  W(B), 
tr(A)  *  tr(B). 


(29) 


Consider  the  matrix 
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A  - 

i  2  i 

. 

°  i  1 

The  visualization  algorithm  produces  the  following  image  for  W(A). 

w  63  sveonds, 
,w(  A)  1 .50000/^. 


Figure  1 


The  image  of  W(A)  in  Rgure  1  is  scaled  by  2,  which  means  that  every  entry  in  A 
is  divided  by  2  before  W(A)  is  computed. 


Next  consider  the  matrix 


0  Jz  ; 


J2i  0  1 


Figure  2 

The  image  of  W(B)  in  Figure  2  is  scaled  by  2. 

The  numerical  ranges  of  A  and  B  are  identical.  We  know  that  W(A)  =  W(B)  is  a 
necessary  but  not  sufficient  condition  for  two  matrices  to  be  unitarily  similar. 
The  matrices  A  and  B  are  not  unitarily  similar.  To  explain  how  these  matrices 
can  have  the  same  numerical  range  and  not  be  unitarily  similar  we  present  the 
following  discussion. 

First  consider  matrix  A.  Define  the  matrix 


’  3 

i 

O' 

'  2 

0 

o' 

'  1  i  O' 

i 

2 

i 

- 

0 

2 

0 

s 

i  0  i 

0 

i 

1 

0 

0 

2 

a 

0  i  -1 

The  matrix  C  has  a  unique  decomposition  into  real  and  imaginary  parts,  Hc  and 
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Kr>  both  symmetric: 


C*Hc  +  iKc 


C  +  C*  ,  C  -  C* 

- — -  4-  j - — - 

2  2i 


1  0  0 


0  1  0 


0  0  0  +1  1  0  1 


0  0  -1 


0  1  0 


For  a  general  matrix,  H  +  iK,  Algorithm  2  graphs  the  support  lines 

x  cos  0  -y  sin  0  a  X(0) 

where  X(0)  is  the  largest  eigenvalue  of  H(6)  »  cos  6  H  -  sin  0  K.  We  compute 
from  above  that  for  C, 


cos  0  Hq-  sin  0  Kq 


cos  0 

-  sin  0 

0 

-  sin  0 

0 

-  sin  0 

0 

-  sin  0 

-COS0 

The  characteristic  polynomial  of  the  preceding  matrix  H  (0)  is 


X3  +  [-2  sin20  -  cos20]X  >  [  -  sin20  cos  0  +  sin20  cos  0]  »  0, 


X3  -  {1  +  sin20]X  »  0. 


Solving  for  X  we  obtain 


X  *  ±  v  1  +  sin  0 
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Thus 


X(0)  s*  -J 1  +  sin20 . 

Now  consider  the  matrix  B.  Note  that  B  is  unitarily  (permutation)  similar  to  the 
matrix 


Define  the  matrix 


B’ 


2  0  0 

0  3  72  j 

0  J2i  1 


D  -  B*  -  2I3 


0  0  0 
0  1  J2i 

0  Jz\  -1 


As  above,  D  *  HD  +  iKD,  for  some  unique  hermitian  HQ  and  KQ: 


and 


D  +  D* 
2 


0  0  0 

0  1  0 

0  0-1 


D-D* 

2i 


0  0 
0  1 
1  0 


To  graph  W(B)  using  the  second  visualization  algorithm,  the  lines 
x  cos  0  -  y  sin  0  «  X(0)  are  graphed  where  X(0)  is  the  largest  eigenvalue  of 


,  ■  s  ,-  w  .**  -  V i  • 
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H(0)  =  cos  9  H  -  sin  9  K  .  For  B', 

Hd(0)  -  cos  0  Hd  -  sin  0  Kd 

'00  0 

»  0  cos©  -^2 sin  0 

0  -  J2  sin  0  -  cos  9 

The  characteristic  polynomial  of  HD(0)  is 

X2  +  [  -  cos20  -  2  sin20  ]  =  0, 


or 


or 


2  2  2  2 
X  +  [  -  cos  0  -  sin  0  -  sin  0  ] 


0, 


X2  *  1  +  sin20. 


Solving  for  X  we  have 

*  X  =  ±V  1  +  sin20 . 

► 

Thus 

X(0)  =  J 1  +  sin20 . 

Hence  the  support  lines  are  the  same  for  W(D)  and  W(C).  But  W(A)  =  W(C)  +  2 
and  W(B) «  W(D)  +  2,  and  thus  W(A)  -  W(B).  The  important  thing  to  note  here  is 
that  the  matrices  A  and  B  have  the  property  that 

X(6)  -  Xmax(cos  0  H  -  sin  0  K) 
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is  the  same  for  all  6  for  both  A  and  B.  If  we  denote  the  maximum  eigenvalue  of 
the  hermitian  part  of  e^A  by  X  (8)  (and  similarly  for  X  (0))  then  the  geometric 

A  B 

condition 

Xa(8)  *  XB(0)  forallQe  [0, 2n] 

is  not  equivalent  to  the  algebraic  condition  that  A  is  unitarily  similar  to  B. 

If  A  and  B  were  unitarily  similar  then  U*AU  =  B  would  imply  that 

u*hau  -  HBi 
and 

U*KaU  =  K0. 

Hence  if  A  were  unitarily  similar  to  B  then  C  ■  A  -  2I3  would  be  unitarily  similar  to 
D  =  B'  -  2I3.  Thus  we  can  work  with  the  matrices  C  and  D  in  our  discussion. 
Consider 

U*HdU  -  Hc  (30) 

and 

U*KdU  =  Kc-  (31) 

From  (30)  we  can  solve  the  system  for  the  matrix  U: 

HdU  =  UHC, 


66 


'0 

0 

0  ' 

'  1 

0 

0  ' 

0 

1 

0 

u  = 

u 

0 

0 

0 

0 

0 

-1 

0 

0 

-1 

- 1 

0 

0 

0 

c 

-A 
— A 

0 

CO 

1 

U21 

u22 

U23 

= 

U21 

0 

■  U23 

/U31 

*  U32 

■  U33 

U31 

0 

*  U33 

The  above  equalities  lead  to  un  »  u13  *  u^  »  u23  =  u32  ■  u31  *  0.  So,  if  a  unitary 
U  exists  satisfying  (30)  and  (31)  then  it  must  have  the  form 


U  = 


U12  0 
0  0 
0  Ugg 


(32) 


From  (31 )  we  have  the  equalities 


KdU  =  UKC, 


0  0  0 
0  0  1 
0  1  0 


0 


0 


U12  0 

0  0 

0  U33 


0  u12  0 

u21  0  0 

0  0  u33 


0 

0 

0 

U12  *  * 

/2 

0 

0 

U33 

= 

0  *  * 

U21 

0 

0 

0  *  * 

0  1  0 
1  0  1 
0  1  0 


After  we  have  computed  the  first  column  of  the  righthand  side  of  the  last 
equality,  we  need  go  no  further.  The  equality  shows  that  u12  ■  0.  This  fact 
combined  with  (32)  contradicts  the  unitary  property  of  U. 
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Open  Questions 

There  is  a  nearly  unlimited  number  of  open  questions  in  this  field.  However, 
much  current  research  is  along  the  following  general  lines. 

Let  A:V  V  be  an  operator  on  a  unitary  space  V.  Let  X  be  a  subset  of  V  and 
/: C  -» C  be  a  complex  function.  Let  M  be  a  subset  of  C.  Describe  the  set 

W(A,  X,  M,  f)  m  {  z  |  /((Ax,  x))  s  M  for  all  x  e  X  }. 


There  are  many  variations  on  this  question.  For  example,  suppose 
W(A,  X,  M,  /)  has  certain  geometric  properties,  e.g.,  symmetry  with  respect  to 
a  line  or  a  point,  then  what  can  be  concluded  about  A? 

Here  are  some  simple  instances  of  this  type  of  question: 


1 .  If  A  is  nilpotent  and  W(A)  is  a  disk,  must  it  be  centered  at  the  origin? 

2.  If  W(A)  is  a  disk  centered  at  the  origin,  is  A  nilpotent?  The  answer  is  "no": 


0  1  0 

0  0  0 

0  0  0.1 


3.  If  A  and  B  are  3  x  3  and  W(A)  *  W(B),  W(A-1)  =  W(B_1)  and  tr(A)  =  tr(B),  is  it 
true  that  A  and  B  are  unitarily  similar?  (W(A)  =  W(B)  is  not  enough  to 
conclude  that  A  and  B  are  unitarily  similar.) 


A  brief  perusal  of  the  bibliography  in  Section  IV  indicates  the  broad  scope  of 
this  research. 
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IIL  Research  of  M.  Marcus,  1 983  - 1 988 


1  .  On  the  equality  of  decomposable  symmetrized  tensors  (with  J.  Chollet), 

Linear  and  Multilinear  Algebra,  13  (1983),  253-266. 

This  paper  continues  earner  work  by  the  authors  on  finding  necessary  and  sufficient 
conditions  for  two  decomposable  symmetrized  tensors  to  be  equal.  In  the  previous  paper 
[Linear  and  Multilinear  Algebra  6  (1978),  317-326]  the  linear  independence  of  the  vectors 
forming  such  tensors  was  assumed.  In  the  present  paper,  this  assumption  is  dropped  and 
much  simpler  requirements  for  equality  are  obtained.  The  paper  also  includes  conditions 
for  a  decomposable  symmetrized  tensor  to  be  0.  This  research  is  related  to  recent  work  of 
JA  Dias  da  Silva,  R.  Merris,  S.  Pierce,  G.N.  de  Oliveira,  and  S.Q.  Williamson. 

2  .  Solution  to  problem  6366,  American  Mathematical  Monthly,  90  (1983, 

409-410. 

3.  Products  of  doubly  stochastic  matrices  (with  K.  Kidman  and  M.  Sandy), 
Linear  and  Multilinear  Algebra,  15  (1984),  331-340. 

In  studying  Westwick's  theorem  on  higher  numerical  ranges,  the  theory  of  elementary 
doubly  stochastic  matrices  arises.  This  concept  is  related  to  the  work  of  M.  Goldberg  and 
E.G.  Straus  [Linear  Algebra  Appl,  18  (1977),  1-24]  on  the  representation  of  a  doubly 
stochastic  matrix  as  a  product  of  elementaiy  doubly  stochastic  matrices.  This  paper  studies 
the  class  of  doubly  stochastic  matrices  that  can  be  written  as  products  of  elementary 
doubly  stochastic  matrices.  The  same  questions  for  orthostochastic  matrices  are  also 
investigated. 

4.  Unitarily  invariant  generalized  matrix  norms  and  Hadamard  product  (with 
K.  Kidman  and  M.  Sandy),  Linear  and  Multilinear  Algebra,  16  (1984),  197- 
213. 

Let  ||  •  ||  be  a  unitarily  invariant  generalized  matrix  norm  on  Mn(C),  the  space  of  n-square 
complex  matrices.  Theorems  are  developed  relating  the  Hadamard  product  (entrywise 
product)  of  two  matrices  A,  B  e  Mn(C)  to  the  singular  values  of  A  and  B.  For  p  2 1 , 

1  s  k  s  n,  let 


where  a  (A)  z  -  z  a  (A)  are  the  singular  values  of  A.  In  this  paper  the  following  inequality 
i  n 


69 


is  proved:  ||  A-B  |£s||A|£  1|B|£.  If  1  <ks  nit  is  also  proved  that  ||  A-B|f-|A|f  ||B|£l 
and  only  if  A  -  a^E,.,  and  B  -  b^E,.,  where  E..  is  the  matrix  with  1  in  position  (i,j)  and  zeros 
elsewhere.  The  case  k- 1  is  also  discussed. 

5.  An  exponential  group,  Linear  and  Multilinear  Algebra,  14  (1984),  293*296. 

Let  M(r)  be  the  n-square  matrix  whose  (i,j)  entry  is 


It  is  proved  that  the  mapping  r  -» M(r)  establishes  an  isomorphism  from  the  additive  group  of 
the  real  numbers  into  the  multiplicative  structure  of  the  n  x  n  matrices.  This  paper 
investigates  M(r)  as  an  exponential  matrix. 

6.  Solution  to  problem  6430  (with  J.  Bruno),  American  Mathematical  Monthly, 
92  (1985),  148-149. 

7.  Conditions  for  the  generalized  numerical  range  to  be  real  (with  M.  Sandy), 
Linear  Algebra  and  Appl.,  71  (1985),  219-239. 

If  A  and  C  are  n-square  complex  matrices  then  the  C-numerical  range  of  A  is  the  totality  of 

numbers  tr(CU*AU)  as  U  varies  over  all  unitaiy  matrices.  This  paper  obtains  necessary  and 
sufficient  conditions  for  the  C-numericaf  range  of  A  to  be  a  subset  of  the  real  axis.  The 
principal  condition  is  that  both  A  and  C  must  be  translates  of  Hermitian  matrices. 

8.  Ryser's  permanent  identity  in  symmetric  algebra  (with  M.  Sandy),  Linear 
and  Multilinear  Algebra,  18  (1985),  183-196. 

The  polynomial  algebra  over  a  field  is  canonically  isomorphic  to  the  symmetric  algebra  over 
a  vector  space.  Several  identities  expressing  homogeneous  polynomials  in  terms  of  sums 
of  powers  of  linear  polynomials  are  exploited  to  obtain  Ryser's  permanent  identity 
[Combinatorial  Mathematics,  MAA  Carus  Monograph  No.  14,  Wiley,  New  York,  1963]  as 
well  as  extensions  of  identities  due  to  Bebiano  [Pacific  J.  Math.,  101  No.  1,  (1982),  1-9] 

9.  Singular  values  and  numerical  radii  (with  M.  Sandy),  Linear  and 
Multilinear  Algebra,  18,  No.  3,  (1985),  337-353. 

The  purpose  of  this  paper  is  to  prove  the  following  result  relating  the  singular  values  and 
the  numerical  radius  of  a  matrix:  For  any  n-square,  complex  matrix  A  with  singular  values 
a  i-io  2  0  and  numerical  radius  r(A) 

in 
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£  «tA). 


ctl  +  -  +  °h 

n 

with  equality  if  and  only  if  A/r(A)  Is  unltarily  similar  to  the  direct  sum  of  a  diagonal  unitary 
matrix  and  unit  multiples  of  2  x  2  matrices  of  the  form 

‘  1  d  ' 

-d  -1 


whereOc  |d|£l. 

10.  Construction  of  orthonormal  bases  in  higher  symmetry  classes  of  tensors 
(with  J.  Chollet),  Linear  and  Multilinear  Algebra,  19  (1986),  133-140. 

A  method  is  presented  for  constructing  an  orthonormal  basis  for  a  symmetry  class  of 
tensors  from  an  orthonormal  basis  of  the  undertying  vector  spaces.  The  basis  so  obtained 
is  not  composed  of  decomposable  symmetrized  tensors.  Indeed,  we  show  that,  for 
symmetry  classes  of  tensors  whose  associated  character  has  degree  higher  than  1 ,  it  is 
impossible  to  construct  an  orthogonal  basis  of  decomposable  symmetrized  tensors  from 
any  basis  of  the  underlying  vector  space.  The  paper  poses  an  open  problem  on  the 
possfcility  of  a  symmetry  class  having  an  orthonormal  basis  of  decomposable  symmetrized 
tensors. 

1 1 .  Computer  generated  numerical  ranges  and  some  resulting  theorems  (with 
C.  Pesce),  Linear  and  Multilinear  Algebra,  21  (1987),  121-157. 

The  numerical  range,  W(A),  of  an  arbitrary  n-square  matrix  A  is  the  union  of  the  numerical 
ranges  of  all  2-square  real  compressions  of  A.  As  a  result,  a  simple  graphics  program  is 
written  that  accurately  exhfoits  W(A)  for  real  A,  and  suggests  several  conjectures  relating 
the  geometry  of  W(A)  to  algebraic  properties  of  A.  Some  of  these  conjectures  are 
analyzed  in  the  final  sections  of  the  paper. 

12.  Solution  to  problem  1231,  Mathematics  Magazine,  60  No.  1  (1987),  42. 

13.  Vertex  points  in  the  numerical  range  of  a  derivation  (with  M.  Sandy), 
Linear  and  Multilinear  Algebra,  21  (1987),  385-394. 

This  paper  contains  a  number  of  results  on  the  distribution  of  values  of  subdeterminants  of 
normal  matrices.  It  is  a  continuation  of  earlier  work  of  M.  Marcus  [Indiana  University  Math.  J., 
22  (1973),  1137-1149]. 
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14.  Solution  to  problem  E3179,  submitted  to  American  Mathematical  Monthly. 

15.  Solution  to  problem  1248,  A  Curious  property  of  1/7,  (with  C.  Pesce), 
Mathematics  Magazine,  60  (1987),  42. 

16.  Two  Determinant  Condensation  Formulas,  Linear  and  Multilinear  Algebra, 
22  (1987),  95-102. 

This  paper  corrects  and  extends  several  classical  results  that  express  the  determinant  of  a 
block  matrix  in  terms  of  determinants  of  the  constituent  blocks. 

17.  Symmetry  properties  of  higher  numerical  ranges  (with  M.  Sandy),  Linear 
Algebra  and  Appl.,  1 04  (1 988),  1 41  -1 64. 

Let  A  be  a  linear  operator  on  a  finite  dimensional  unitary  space  V  of  dimension  n.  The  kth 
higher  numerical  range  of  A,  denoted  by  Wk(A),  is  the  totality  of  complex  numbers  tr(PAP) 
where  P  runs  over  ail  k-dimensional  orthogonal  projections  on  V.  It  is  proved  that  Wk(A)  is 

polygon  with  the  real  axis  as  a  line  of  symmetry,  k  - 1 . n,  if  and  only  if  A  is  normal  with  a 

real  characteristic  polynomial.  Several  non-normal  examples  are  exhibited  that  reveal  the 
extent  to  which  the  symmetry  of  all  the  Wk(A)  is  required. 

18.  Advanced  problem,  Triangular  Kronecker  Products,  (with  C.  Pesce), 
accepted  for  publication,  American  Math.  Monthly. 

19.  Bessel's  inequality  in  tensor  space  (with  M.  Sandy),  Linear  and  Multilinear 
Algebra,  in  press. 

Let  A  be  an  n-square  complex  matrix  and  define  AA  to  be  the  nl-square  matrix  whose 
entries  are 


n 

i-t 


a 


o(i),T(i) 


where  a  and  t  run  lexicographically  over  Sn>  If  A  is  positive  definite  Hermitian  and  %  is  a  unit 
nl-tuple  then 


(AAx-x)2det(A)  +  lXx(oM2c(A) 
where  c(A)  is  the  largest  of  the  numbers 
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and  the  summation  is  over  a  e  S  .  For  n  -  3,  if  A  is  not  permutation  similar  to  a  direct  sum 

n 

and  x  is  a  unit  n!-tuple  then  (A^%,  x)  ■  det(A)  iff  x  is  a  multiple  of  the  alternating  character. 
The  relationships  among  recent  results  of  Bapat  and  Sunder,  Chollet,  and  Gregorac  and 
Henzel  are  also  discussed. 

20.  A  unified  approach  to  some  classical  matrix  theorems,  submitted. 

An  elementary  inequality  is  proved  that  obtains  the  lower  bound  of  the  product  of  forms 

(Axx)  (A*1x,x),  where  x  is  a  unit  vector  and  A  is  a  positive  definite  Hermitian  matrix.  Using 
this  inequality  it  is  possible  to  provide  a  unified  treatment  of  the  following  theorems:  the 
Hadamard  determinant  theorem;  the  Fischer  inequality;  the  Kantorovich  inequality;  Weyl's 
inequalities. 

Books 

1.  Discrete  Mathematics:  A  Computational  Approach  Using  BASIC, 
Computer  Science  Press,  Rockville,  Maryland,  1983. 

2.  MacAlgebra,  Basic  Algebra  on  the  Macintosh  (with  R.  Marcus  and  C. 
Baczynski),  Computer  Science  Press,  Rockville,  Maryland,  1986. 

3.  An  Introduction  to  Pascal  and  Precalculus,  Computer  Science  Press, 
Rockville,  Maryland,  1986. 

4.  Computing  Without  Mathematics:  BASIC,  Pascal,  Applications  (with  J. 
Marcus),  Computer  Science  Press,  Rockville,  Maryland,  1986. 

5.  Introduction  to  Linear  Algebra  (with  H.  Mine),  Dover  Publications  Inc.,  New 
York,  1988. 

6.  Foundations  of  Numerical  Linear  Algebra  I,  Center  for  Computational 
Sciences  and  Engineering,  University  of  California,  Santa  Barbara, 
Monographs  in  Scientific  Computation,  1988. 


Science  Citation  Index 


A  search  was  recently  conducted  of  the  Science  Citation  Index  for  the  total 
number  of  references  to  the  work  of  M.  Marcus  since  1983.  Self  references 
were  excluded  in  the  search  criteria.  The  number  of  such  citations  is  593. 


IV.  Numerical  Range  Bibliography 

The  following  pages  contain  a  printout  of  the  Numerical  Range  Bibliography.  It 
is  sorted  alphabetically  by  the  first  author. 
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V.  Appendix 


VITA 

PERSONAL  BACKGROUND 

Born:  Albuquerque,  New  Mexico,  July  31 ,  1927 

Education:  Attended  public  schools  in  California 

Military  Service:  United  States  Navy,  1944-1946, 

honorable  discharge 

Married:  Rebecca  Elizabeth  Marcus 

Children:  Jeffrey  (employed,  Micropoint,  Los  Angeles) 

Karen  (Ph.D.  student,  Stanford  University) 

Academic  Degrees: 

1950  B.A.  (highest  honors  in  Mathematics)  University  of 

California,  Berkeley 

1 953  Ph.D,  University  of  California,  Berkeley 

PROFESSIONAL  EXPERIENCE 

1 987  -  present  Professor  of  Computer  Science,  UCSB 

1983-  1987  Professor  of  Mathematics  and  Computer  Science, 

University  of  California,  Santa  Barbara 

1979 -present  Founder,  Microcomputer  Laboratory, University  of 

California,  Santa  Barbara 

1978-  1986  Associate  Vice  Chancellor  and  Dean,  Research  and 
Academic  Development,  University  of  California,  Santa 
Barbara 
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1973  -  1979  Director,  Institute  for  the  Interdisciplinary  Applications  of 
Algebra  and  Combinatorics, University  of  California,  Santa 
Barbara 

1963-  1968  Chairman,  Department  of  Mathematics,  University  of 

California,  Santa  Barbara 

1962  -  1983  Professor  of  Mathematics,  University  of  California,  Santa 
Barbara 

1960-  1961  Research  Mathematician,  U.S.  National  Bureau  of 

Standards,  Washington,  D.C. 

1954-  1961  Instructor,  Assistant  and  Associate  Professor  of 

Mathematics,  University  of  British  Columbia 

UNDERGRADUATE  ACADEMIC  ACTIVITIES 

1964-  1972  Lecturer,  National  Science  Foundation  Linear  Algebra 

Conference  for  College  Teachers,  University  of  California, 
Santa  Barbara 

1965  -  1966  Lecturer,  National  Science  Foundation  In-service  Institute 
for  Secondary  School  Teachers,  University  of  California, 
Santa  Barbara 

1965-  1974  Visiting  lecturer  for  the  Mathematical  Association  of 

America,  touring  four-year  undergraduate  institutions  in  the 
far  western  area  giving  lectures  on  undergraduate 
mathematics 

1965  -  present  Author  and  co-author  of  20  undergraduate  textbooks  (see 
publication  list) 

1975  Principal  Investigator,  Summer  Projects  Grant  and 

Regents'  Undergraduate  Instructional  Improvement  Grant 
for  training  Scientific  Information  Specialists,  University  of 
California,  Santa  Barbara 
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1979  -  1986 


Established  Microcomputer  Laboratory  at  the  University  of 
California  at  Santa  Barbara,  under  grants  from:  the  Fund 
for  the  Improvement  of  Postsecondary  Education; 
California  Postsecondary  Education  Commission; 
Instructional  Scientific  Equipment  Program,  NSF. 

1979  -  1984  Principal  Investigator,  The  Comprehensive  Program,  Fund 
for  the  Improvement  of  Postsecondary  Education, 
Curriculum  Development  Project  in  Applied  Algebra 

1979-  1984  Program  Director,  Intensive  Short  Course  in  Basic  College 
Level  Mathematics  for  Adult  Reentry  Women  under  grants 
from  the  California  Postsecondary  Education  Commission 
and  The  Development  in  Science  Education  Project  of  the 
National  Science  Foundation 

1987  -  1988  Principal  Investigator,  National  Science  Foundation  Grant, 
Computing  and  Algorithmic  Mathematics  for  Secondary 
School  Teachers 


GRADUATE  ACADEMIC  ACTIVITIES 

1 964  -  present  Author  and  co-author  of  four  graduate  textbooks 

1970  Ford  Foundation  Visiting  Distinguished  Professor, 
University  of  Islamabad,  Islamabad,  Pakistan;  Consultant 
on  curriculum  design  at  the  new  Pakistan  National 
University 

1971  Visiting  Lecturer,  University  of  Victoria,  Victoria,  British 
Columbia;  Assist  in  the  graduate  program 

1973,  1977  Director,  Conferences  on  Matrix  Theory,  sponsored  by  the 
National  Science  Foundation,  University  of  California, 
Santa  Barbara 

1974  Visiting  Distinguished  Professor,  Laval  University, 

Quebec,  Canada;  assist  in  the  graduate  program 
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GRADUATE  STUDENTS 


The  following  mathematicians  have  completed  their  Ph.D.  work  under 
the  direction  of  M.  Marcus: 

Dr.  Roy  Westwick,  Professor 
University  of  British  Columbia 
Vancouver,  B.C.,  Canada 

Thesis:  Linear  transformations  of  Grassmann  algebras 

1960 

Dr.  Nisar  A.  Khan,  Professor 
Muslim  University,  Aligarh,  India 
Thesis:  Matrix  commutators 

1961 

Dr.  Peter  Botta,  Assoc.  Professor 

University  of  Toronto 

Toronto,  Ontario,  Canada 

Thesis:  Linear  transformations  on  algebras 

1965 

Dr.  Stanley  G.  Williamson,  Professor 

University  of  California 

San  Diego,  California 

Thesis:  Tensor  Algebras 

1965 

Dr.  William  R.  Gordon,  Professor 
Department  of  Mathematics 
University  of  Victoria 
Victoria,  B.C.,  Canada 

Thesis:  Inequalities  for  generalized  matrix  functions 

1965 

Dr.  George  Soules 
Institute  for  Defense  Analysis 
Princeton,  New  Jersey 
Thesis:  Combinatorial  functions 

1966 
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Dr.  Paul  J.  Nikolai,  Mathematician 
Wright-Patterson  Air  Force  Base 

Thesis:  Mean  value  properties  of  generalized  matrix  functions 

(This  thesis  was  supervised  jointly  with  Professor  H.  J.  Ryser  (deceased), 

California  Institute  of  Technology) 

1966 

Dr.  Stephen  J.  Pierce,  Professor 
California  State  University 
San  Diego 

Thesis:  Generalized  isometries 

1968 

Dr.  William  Watkins,  Professor 
California  State  University,  Northridge 
Northridge,  California 

Thesis:  Inequalities  for  derivation  operators  on  a  tensor  space 

1969 

Dr.  Russell  Merris,  Professor 
California  State  University  at  Hayward 
Hayward,  California 

Thesis:  A  generalization  of  the  associated  transformation 

1969 

Dr.  Mohammad  Shafqat  Ali,  Assoc.  Professor 
California  State  University  at  Long  Beach 
Long  Beach,  California 

Thesis:  Additive  commutators,  Jordan  products  and  bilinear  functions 

1970 

Dr.  Elizabeth  Wilson,  Mathematician 

Naval  Labs.  Pt.  Mugu,  California 

Thesis:  Partial  derivations  on  symmetry  classes  of  tensors 

1971 
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Dr.  James  Holmes,  Assistant  Professor 
Westmont  College 
Santa  Barbara,  California 

Thesis:  Application  of  derivations  to  invariance  problems 
1971 

Dr.  Herbert  Robinson,  Professor 
Department  of  Mathematics 
Texas  A  &  M  University 
College  Station,  Texas 

Thesis:  Quadratic  &  bilinear  forms  on  symmetry  classes  of  tensors 

1975 

Dr.  Patricia  Andresen 
University  of  Alaska 
Fairbanks,  Alaska 

Thesis:  The  finite  dimensional  numerical  range 

1976 

Dr.  Robert  Grone 
University  of  Auburn 
Auburn,  Alabama 

Thesis:  Isometries  of  Matrix  Algebras 

1976 

Dr.  Ivan  Filippenko,  Research  Mathematician 
Lockheed  Aircraft 
Los  Angeles,  California 

Thesis:  Higher  and  Decomposable  Numerical  Ranges 

1977 

Dr.  John  Chollet,  Assistant  Professor 

University  of  British  Columbia 

Vancouver,  British  Columbia,  Canada 

Thesis:  Equalities  of  decomposable  symmetrized  tensors 

1979 


Dr.  Kenneth  Moore 
Radar  Systems  Group 
Hughes  Aircraft  Co. 

El  Segundo,  California 

Thesis:  Determinantal  Inequalities 

1980 

Dr.  Kent  Kidman 
Hughes  Aircraft  Co. 

El  Segundo,  California 

Thesis:  Stochastic  Matrices  and  unitarily  Invariant  Norms 
1983 

Claire  Pesce 

Naval  Weapons  Center,  China  Lake 
China  Lake,  California 

Thesis:  Visualization  of  the  Numerical  Range 
1988 


ACADEMIC  AWARDS  AND  DISTINCTIONS 

1950  Graduated  highest  honors  in  mathematics,  University  of 

California,  Berkeley 

1 954  Fulbright  Award 

1956-57  National  Research  Council,  National  Science  Foundation, 

Post-doctoral  Research  Fellowship 

1956,  1958-60,  National  Science  Foundation  Research  Grants 
1975-84 

1962  Certificate  of  Award  for  Distinguished  Service,  U.S. 

Department  of  Commerce,  National  Bureau  of  Standards 

1962 -present  Principal  Investigator  on  Air  Force  Office  of  Scientific 
Research  grants 
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1 965  Mathematical  Association  of  America  Editorial  Prize  for  the 
article  entitled:  "Linear  Transformations  on  matrices" 

1966  L.R.  Ford  Memorial  Prize  awarded  by  the  Mathematical 
Association  of  America  for  the  article,  "Permanents" 


EDITORIAL  ACTIVITIES 

1.  Mathematics  Editor,  Computer  Science  Press 

2.  Editor,  Linear  and  Multilinear  Algebra,  published  by  Gordon  and  Breach, 
Science  Publishers  Inc. 

3.  Associate  Editor,  Linear  Algebra  and  Its  Applications,  Elsevier  Science 
Publishing  Co.,  Inc. 

4.  Member  of  the  Editorial  Board,  Pure  and  Applied  Mathematics  Series, 
Marcel  Dekker,  Inc. 

5.  Editor,  Linear  Algebra  Volumes  of  Encyclopedia  of  Applicable 
Mathematics,  Addison-Wesley  Publishing  Co. 

6.  Member  of  the  Editorial  Board,  Linear  Algebra  and  Its  Applications 

7.  Associate  Editor,  Advanced  Problem  Section,  American  Mathematical 
Monthly 

8.  Referee  and  Reviewer  for  the  following  journals: 

Linear  and  Multilinear  Algebra 
Linear  Algebra  and  Its  Applications 
Duke  Journal 
Proceedings  of  the  AMS 
Transactions  of  the  AMS 
Bulletin  of  the  AMS 
Mathematical  Reviews 
Memoirs  of  the  MAA 
American  Mathematical  Monthly 
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Canadian  Journal  of  Mathematics 

Pacific  Journal  of  Mathematics 

Proceedings  of  the  Cambridge  Philosophical  Society 

Zentralblatt 

9.  Technical  reviewer  for  the  Air  Force  Office  of  Scientific  Research 

10.  Technical  reviewer  for  the  Mathematics  Division  of  the  National  Science 
Foundation 

1 1 .  Technical  reviewer  for  the  National  Research  Council  of  Canada 

12.  Technical  reviewer  for  United  States-lsrael  Binational  Science 
Foundation 

1 3.  Editorial  advisor  for  the  following  publishers: 

Houghton-Mifflin  Company 

W.A.  Benjamin,  Inc. 

Harcourt,  Brace  and  World 

1 4.  Advisory  Editor,  Letters  in  Linear  Algebra 

15.  Editorial  Board,  Algebras,  Groups,  and  Geometries 


SELECTED  INVITED  PAPERS 

1963  International  Conference,  "Recent  Advances  in  Matrix  Theory",  U.S. 
Army  Research  Center,  Madison,  Wisconsin 

1 965  Far-Western  meeting  of  the  Mathematical  Association  of  America 

1965  Invited  speaker,  Annual  meeting  of  the  American  Mathematical 
Society 

1965,  1967,  1969 

Symposium  on  Inequalities,  sposored  by  Aerospace  Resarch 
Laboratories,  U.S.  Air  Force 
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1967  International  Symposium  on  Combinatorial  Analysis,  sponsored  by 
the  Society  for  Industrial  and  Applied  Mathematics 

1 972  Conference  on  Numerical  Algebra,  Los  Alamos  Scientific  Laboratory 

1 974  University  of  California,  Los  Angeles 

1975  University  of  Chicago,  Chicago,  Illinois 

1 975  University  of  California,  San  Diego 

1 975  California  State  University  of  Hayward 

1975  AMS  meeting,  Kalamazoo,  Michigan,  Special  Session  on  Matrix 
Theory 

1 975  California  Mathematical  Council  Conference,  Asilomar,  California 

1976  Northern  California  Section  of  the  MAA  annual  meeting.  University  of 
California,  Davis 

1977  Gatlinburg  VII,  Conference  on  Numerical  Algebra 

1978  New  York  Academy  of  Science,  Second  International  Conference  on 
Combinatorial  Mathematics 

1980  California  Institute  of  Technology  Colloquium  Series 

1 980  Oberwolfach  Conference  on  General  Inequalities 

1981  Conference  on  Numerical  Algebra,  Oxford  University 

1982  Mid  Atlantic  Conference  on  Educational  Computing,  Bennett  College, 
Greensboro,  N.C. 

1984  Invited  contribution  Special  Issue  of  Linear  Algebra  and  Its 
Applications  honoring  Helmut  Wielandt 

1 986  University  of  California,  Riverside 
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1 986  University  of  California,  San  Diego 

1 986  Conference  on  Computers  and  Mathematics,  Stanford  University 

1 986  Western  Educational  Computing  Consortium,  Irvine 

1988  SIAM  Conference  on  Applied  Linear  Algebra,  Madison  Wisconsin 

MEMBERSHIP  IN  LEARNED  SOCIETIES 

American  Mathematical  Society 
Mathematical  Association  of  America 
American  Association  of  University  Professors 
Sigma  Psi;  Pi  Mu  Epsilon 

American  Association  for  the  Advancement  of  Science 

Washington  Academy  of  Science 

Society  for  Industrial  and  Applied  Mathematics 

Society  for  Technical  Communication 

Association  for  Computing  Machinery 

UCSB  UNIVERSITY  SERVICE 

Associate  Vice  Chancellor,  Research  and  Academic  Development.  (1979  - 
1987) 

The  following  units  reported  to  this  office: 

ACTER 

Intructional  Development/Leaming  Resources 
Microcomputer  Laboratory 
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Off  Campus  Studies 

University  Center  at  Ventura 

University  Extension 

Algebra  Institute 

Center  for  Black  Studies 

Center  for  Chicano  Studies 

Community  and  Organization  Research  Institute 

Computer  Systems  Laboratory 

Intercampus  Institute  for  Research  at  Particle  Accelerators 

Institute  for  Polymers  and  Organic  Solids 

Institute  of  Environmental  Stress 

Marine  Science  Institute 

Quantum  Institute 

Social  Process  Research  Institute 


Numerous  ad  hoc 
1962  - 1963 
1962  - 1963 

1962-  1963 
1964-  1965 

1963-  1964 

1969-  1972 
1975-  1977 

1970  -  1972 

1973-  1974 

1973-  1974 
1974 

1974-  1975 


personnel  review  committees 
Chairman,  Statistics  Committee 
Academic  Senate  Educational  Policy  Committee 
Chairman,  Computer  Committee 

Digital  Computer  Committee 
Academic  Senate  Research  Committee 

Academic  Senate  Education  Abroad  Committee 
Chairman,  Undergraduate  Committee 
Computer  Science  Laboratory  Director  Search  Committee 
Chancellor's  task  force  on  career  development 
Academic  Senate  Athletic  Policy  Committee 
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1975 


Ad  Hoc  Committee  for  Scientific  Communication 


1979-  1986 


Coordinator,  Chinese  Exchange  Program 


1988  -  Present  Computing  Task  Force 


PROFESSIONAL  REFERENCES 

Professor  Stephen  P.  Diliberto 
Department  of  Mathematics 
University  of  California 
Berkeley,  California  94720 
(415)  642-6550 

Professor  Ky  Fan,  Emeritus 
Department  of  Mathematics 
University  of  California 
Santa  Barbara,  California  93106 
(805)  961-2171 

Professor  Marshall  Hall 
Department  of  Mathematics 
Emory  University 
Atlanta,  Georgia  30322 
(404)  727-5605 

Raymond  Huerta,  Coordinator 
Affirmative  Action 
University  of  California 
Santa  Barbara,  California  93106 
(805)  961-2089 

Professor  Robert  A.  Huttenback 
Former  Chancellor 
2661  Todos  Santos 
Santa  Barbara,  California  93105 
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Professor  Robert  Mehrabian 
Dean,  College  of  Engineering 
University  of  California 
Santa  Barbara,  California  93106 
(805)  961-3141 

Professor  William  H.  Meyer,  Emeritus 
Administrative  Head 
Department  of  Mathematics 
University  of  Chicago 
Chicago,  Illinois  60637 
(312)  962-7100 

Professor  Henryk  Mine 
Department  of  Mathematics 
University  of  California 
Santa  Barbara,  California  93106 
(805)  961-2171 

Professor  Benjamin  N.  Moyls,  Emeritus 
Department  of  Mathematics 
University  of  British  Columbia 
Vancouver,  British  Columbia,  CANADA 
(604)  228-2848 

Professor  Morris  Newman 
Department  of  Mathematics 
University  of  California 
Santa  Barbara,  California  93106 
(805)  961-2171 

Professor  Gian-Carlo  Rota 
Department  of  Mathematics 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 
(617)  253-4381 
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Professor  Hans  Schneider 
Department  of  Mathematics 
University  of  Wisconsin 
Madison,  Wisconsin  53706 
(608)  263-3053 

Professor  David  S.  Simonett 
Dean,  Graduate  Division 
University  of  California 
Santa  Barbara,  California  93106 
(805)  961-2013 

Professor  David  A.  Sprecher 
Dean,  College  of  Letters  and  Science 
University  of  California 
Santa  Barbara,  California  93106 
(805)  961-3506 

Professor  Olga  Taussky-Todd 
Department  of  Mathematics 
California  Institute  of  Technology 
Pasadena,  California  91125 
(818)  356-4332 

Professor  Stanley  G.  Williamson 
Department  of  Mathematics 
University  of  California 
at  San  Diego 

La  Jolla,  California  92037 
(714)  452-3590 


i 


EXTRAMURAL  SUPPORT 
Air  Force 

Multilinear  Algebra 
M.  Marcus,  H.  Mine 
10/01/66  -  09/30/67 

Air  Force 

Multilinear  Algebra 
M.  Marcus,  H.  Mine 
10/01/67  -  09/30/68 

Air  Force 

Multilinear  Algebra 
M.  Marcus,  H.  Mine 
10/01/68  -  09/30/69 

Air  Force 

Multilinear  Methods 
M.  Marcus,  H.  Mine 
10/01/69  -  09/30/70 

Air  Force 

Multilinear  Methods 
M.  Marcus,  H.  Mine 
10/01/70  -  09/30/71 

Air  Force 

Eigenvalue  Investigators  and  Stability 
R.C.  Thompson,  M.  Marcus,  H.  Mine 
10/01/71  -  09/30/72 

Air  Force 

Inequalities,  Combinatorics  and  Applications 
M.Marcus,  H.  Mine,  R.C.  Thompson 
10/01/72  -  09/30/73 


5' 


mlm 


$51,436 


$61,610 


$60,297 


$62,270 


$60,416 


$44,606 


$42,961 


•r'M* 


< 


Air  Force 

The  Algebraic  Eigenvalue  Problem  with  Applications 
M.  Marcus,  H.  Mine,  R.C.  Thompson 
10/01/73  -  09/30/74 

National  Science  Foundation 
Theoretical  Matrix  Theory 
M.  Marcus 
10/15/73-10/14/74 

Air  Force 

Supplementary  Request 
M.  Marcus,  H.  Mine,  R.C.  Thompson 
06/30/74  -  09/30/74 

Air  Force 

Algebraic  Stability:  A  Linear  Algebra  Bibliography 
M.  Marcus,  H.  Mine,  R.C.  Thompson 
10/01/74  -  09/30/75 

National  Science  Foundation 
Undergraduate  Research  Participation 
M.  Marcus 
02/15/75  -  05/31/76 

Air  Force 

Foundations  of  Stability,  Linear  Algebra  Bibliography 
M.  Marcus,  H.  Mine,  R.C.  Thompson 
10/01/75  -  09/30/76 

National  Science  Foundation 
Computer  Searchable  Information  Files 
M.  Marcus 
07/01/76  - 12/31/78 

Air  Force 

Eigenvalue  Problems  in  Stability  Theory 
M.  Marcus,  R.C.  Thompson,  H.  Mine 
10/01/76  -  09/30/77 


$36,079 


$12,100 


$  7,435 


$55,183 


$26,740 


$51 ,702 


$68,348 


$58,793 
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National  Science  Foundation 
Dissemination  of  Scientific  Information 
M.  Marcus 
09/07/77  -  01/31/79 

Air  Force 

Supplement  to:  The  Localization  of  Eigenvalues 
M.  Marcus 
10/01/77  -  09/30/78 

National  Science  Foundation 
Research  Conference  on  Linear  Algebra 
M.  Marcus 
11/01/77-10/31/78 

Air  Force 

Stability,  Control  and  Numerical  Linear  Algebra 

M.  Marcus 

10/01/78-09/30/79 

International  Business  Machines  Corp. 

Intensive  Short  Course  in  Basic  College  Mathematics 
M.  Marcus 
08/01/79  -  09/30/80 
Air  Force 

Foundations  of  Eigenvalue  Distribution  Theory 
M.  Marcus  et  al 
09/30/79  -  09/29/80 

Cal  Post  Secondary  Education  Commission 
Intensive  Short  Course  in  Basic  College  Mathematics 
M.  Marcus 
10/01/79  -  06/30/80 

Cal  Post  Secondary  Education  Commission 
Intensive  Short  Course  in  Basic  College  Level  Math 
M.  Marcus 
07/01/80  -  09/30/81 


$10,000 


$102,881 


$  2,800 


$80,949 


$  5,000 


$84,143 


$33,000 


$40,000 


Air  Force 

Eigenvalue  Localization  Techniques  in  Numerical  Algebra 
M.  Marcus  et  al 

09/30/80  -  09/30/81  $1 01 ,993 

National  Science  Foundation 

Microcomputer  Equipment  for  Undergraduate  Applied  Mathematics 
M.  Marcus 

1 0/1 5/80  -  09/30/83  $  1 9,31 9 

National  Science  Foundation 

Intensive  Computer  Based  Mathematics  Training 

M.  Marcus 

03/01/81  -  10/31/83  $192,012 

National  Science  Foundation 

Research  Conference  on  Multilinear  Algebra 

R.  Merris,  M.  Marcus 

03/15/81  -08/31/81  $7,550 

Department  of  Education 
A  Program  In  Quantitative  Decision  Making 
M.  Marcus 

09/1 5/81  -  09/1 4/84  $113,961 

Air  Force 

Eigenvalues,  Numerical  Ranges,  Stability  Analysis 
M.  Marcus  et  al 

09/30/81  -  04/30/83  $1 1 4,545 

Air  Force 

Questions  in  Numerical  Analysis  /  Associated  Problems 
M.  Marcus,  M.  Goldberg 

05/01/83  -  04/30/84  $57,51 5 

Department  of  Education 
A  Program  In  Quantitative  Decision  Making 
M.  Marcus 

05/01/83-09/14/84  $6,035 


105 


$56,953 


Air  Fores 

Stability  Analysis  of  Finite  Difference  Schemes 
M.  Marcus,  M.  Goldberg 
05/01/84  -  04/30  '85 

Air  Force 

Stability  Analysis  of  Finite  Difference  Schemes 
M.  Marcus,  M.  Goldberg 
05/01/85  -  04/30/86  $64,444 

Air  Force 

Stability  Analysis  of  Finite  Difference  Schemes 
M.  Marcus,  M.  Goldberg 

05/01/86  -  04/30/87  $75,030 

Air  Force 

Stability  Analysis  of  Finite  Difference  Approximations  to  Hyperbolic  Systems, 
and  Problems  in  Applied  and  Computational  Linear  Algebra 
M.  Marcus,  M.  Goldberg 

5/1/8787  -  4/30/88  $73,693 

National  Science  Foundation 

Computing  and  Algorithmic  Mathematics  for  Secondary  School  Teachers 
M.  Marcus,  J.  Bruno 

3/1 1/87  -  8/31/89  $516,999 

National  Science  Foundation 

A  National  Institute  for  Secondary  School  Teachers  for  the  Dissemination  of 
Computer  Science  and  Algorithmic  Mathematics 
M.  Marcus,  J.  Bruno,  R.  Mayer 

2/1/88  -  8/31/89  $509,560 

Air  Force 

Stability  Analysis  of  Finite  Difference  Approximations  to  Hyperbolic  Systems, 
and  Problems  in  Applied  and  Compuational  Matrix  Theory 
M.  Marcus,  M.  Goldberg 

5/1/88-10/31/88  $78,114 

TOTAL  $3,046,472 


Publication  List 


RESEARCH  PAPERS 


1955 

1.  Reid  convexity  of  a  square  matrix  (with  B.N.  Moyls),  Proc.  Amer.  Math. 
Soc.  6  (1955),  981-983. 

2.  A  remark  on  a  norm  inequality  for  square  matrices,  Proc.  Amer.  Math. 
Soc.  6  (1955),  117-119. 

3.  Some  results  on  the  asymptotic  behavior  of  linear  systems,  Canad.  J. 
Math.  7  (1955),  531-538. 

4.  Boundedness  of  a  continuous  function,  Amer.  Math.  Monthly  62  (1955). 

1956 

5.  An  invariant  surface  theorem  for  a  non-degenerate  system. 
Contributions  to  non-linear  oscillations,  Annals  of  Math.  Study  36 
(1956),  243-256. 

6.  A  note  on  the  existence  of  periodic  solutions  of  differential  equations 
(with  S.P.  Diliberto),  Annals  of  Math.  Study  36  (1956),  237-241. 

7.  Repeating  solutions  for  a  degenerate  system,  Annals  of  Math.  Study  36 
(1956),  261-268. 

8.  On  the  optimum  gradient  method  for  systems  of  linear  equations,  Proc. 
Amer.  Math.  Soc.  1  (1956),  77-81. 

9.  Extramural  properties  of  Hermitian  matrices  (with  J.  McGregor),  Canad. 
J.  Math.  8  (1956),  524-531. 

10.  An  eigenvalue  inequality  for  the  product  of  normal  matrices,  Amer. 
Math.  Monthly  63  (1956),  173-174. 
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1957 


11.  On  the  maximum  principle  of  Ky  Fan  (with  B.  Moyls),  Canad.  J.  Math.  9 
(1957),  313-320. 

12.  Inequalities  for  symmetric  functions  and  Hermitian  matrices  (with  L. 
Lopes),  Canad.  J.  Math.  9  (1957),  304-312. 

13.  A  note  on  symmetric  functions  of  eigenvalues  (with  R.C.  Thompson), 
Duke  Math.  J.  24  (1957),  43-46. 

14.  A  note  on  the  values  of  a  quadratic  form,  J.  Wash.  Acad.  Sci.  47  (1 957), 
97-99. 

15.  Some  extreme  value  results  for  indefinite  Hermitian  matrices  I.  (with  B. 
Moyls  and  R.  Westwick),  Illinois  J.  Math.  1  (1957),  449-457. 

16.  On  subdeterminants  of  doubly  stochastic  matrices,  Illinois  J.  Math.  1 
(1957),  583-590. 

17.  A  determinantal  inequality  of  H.P.  Robertson  II,  J.  Wash.  Acad.  Sci.  47 
(1957),  264-266. 

1 8.  Maximum  and  minimum  values  for  the  elementary  symmetric  functions 
of  Hermitian  forms  (with  B.  Moyls),  J.  Lond.  Math.  Soc.  32  (1957),  375- 
377. 

19.  Convex  functions  of  quadratic  forms,  Duke  J.  Math.  24  (1957),  321- 
326. 

1958 

20.  Some  extreme  value  results  for  indefinite  Hermitian  matrices  II,  (with  B. 
Moyls  and  R.  Westwick),  Illinois  J.  Math.  2  (1958),  408-414. 

21.  On  a  determinantal  inequality,  Amer.  Math.  Monthly  65  (1958),  266- 
268. 
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22.  On  doubly  stochastic  transforms  of  a  vector,  Quart.  J.  Math.  Oxford  2 
(1958),  74-80. 

1959 

23.  On  the  minimum  of  the  permanent  of  a  doubly  stochastic  matrix  (with  M. 
Newman),  Duke  J.  Math.  26  (1959),  61-72. 

24.  Convexity  of  the  field  of  a  linear  transformation  (with  A.  Goldman), 
Canad.  Math.  Bull.  2  (1959),  15-18. 

25.  Linear  transformations  on  algebras  of  matrices  (with  B.  Moyls),  Canad. 
J.  Math.  11  (1959),  61-66. 

26.  All  linear  operators  leaving  the  unitary  group  invariant,  Duke  J.  Math. 
26  (1959),  155-163. 

27.  Extremal  properties  of  Hermitian  matrices  II  (with  B.  Moyls  and  R. 
Westwick),  Canad.  J.  Math.  11  (1959),  379-382. 

28.  Linear  transformations  on  algebras  of  matrices  II  (with  R.  Purves), 
Canad.  J.  Math.  11  (1959),  383-396. 

29.  A  note  on  the  Hadamard  product  (with  N.  Khan),  Canad.  Math.  Bull.  2 
(1959),  81-83. 

30.  Transformations  on  tensor  product  spaces  (with  B.  Moyls),  Pacific  J. 
Math.  9(1959),  1215-1221. 

31 .  Diagonals  of  doubly  stochastic  matrices  (with  R.  Ree),  Oxford  Quart.  J. 
Math.  10  (1959),  296-302. 

1960 

32.  On  matrix  commutators  (with  N.  Khan),  Canad.  J.  Math.  12  (1960),  269- 
277. 

33.  Space  of  k-commutative  matrices  (with  N.  Khan),  J.  Research  Nat’l 
Bureau  of  Standards  64B  (1960),  51-54. 
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34.  Some  properties  and  applications  of  doubly  stochastic  matrices,  Amer. 
Math.  Monthly  67  (1960),  215-220. 

35.  A  note  on  a  group  defined  by  a  quadratic  form  (with  N.  Khan),  Canad. 
Math.  Bull.  3  (1960),  143-148. 

36.  Linear  maps  on  skew-symmetric  matrices;  the  invariance  of  elementary 
symmetric  functions  (with  R.  Westwick),  Pacific  J.  Math.  10  (1960),  917- 
924. 

37.  The  maximum  number  of  equal  nonzero  subdeterminants  (with  H. 
Mine),  Archiv.  D.  Math.  11  (1960),  95-100. 

38.  Permanents  of  doubly  stochastic  matrices  (with  M.  Newman),  Proc.  of 
Symposia  in  Applied  Math.  10,  Amer.  Math.  Soc.,  1960. 

39.  On  a  commutator  result  of  Taussky  and  Zassenhaus  (with  N.  Khan), 
Pacific  J.  Math.  10  (1960)  1337-1346. 

40.  On  a  theorem  of  I.  Schur  concerning  matrix  transformations  (with  F. 
May),  Archiv.  D.  Math.  11  (1960),  401-404. 

1961 

41.  On  the  unitary  completion  of  a  matrix  (with  P.  Greiner),  Illinois  J.  Math. 
5(1961)  152-158. 

42.  Some  generalizations  of  Kantorovich's  inequality  (with  N.  Khan), 
Portugal.  Math.  20  (1961),  33-38. 

43.  Another  extension  of  Heinz's  inequality,  J.  of  Research  Nat'l  Bureau  of 
Standards  65B  (1961),  129-130. 

44.  A  note  on  normal  matrices  (with  N.  Khan),  Canad.  Math.  Bull.  4  (1961), 
23-27. 

45.  Symmetric  means  and  matrix  inequalities  (with  P.  Bullen),  Proc.  Amer. 
Math.  Soc.  12  (1961),  285-290. 
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46.  Bound  for  the  p-condition  number  of  matrices  with  positive  roots  (with 
P.  David  and  E.  Haynsworth),  J.  of  Research,  Nat’l  Bureau  of 
Standards  65  (1961),  13-14. 

47.  The  permanent  function  as  an  inner  product  (with  M.  Newman),  Bull. 
Amer.  Math.  Soc.  67  (1961),  223-224. 

48.  Comparison  theorems  for  symmetric  functions  of  characteristic  roots,  J. 
of  Research  Nat'l  Bureau  of  Standards  65  (1961),113-116. 

49.  Some  results  on  non-negative  matrices  (with  H.  Mine  and  B.  Moyls),  J. 
of  Research  NafI  Bureau  of  Standards  65  (1961),  205-209. 

50.  On  the  relation  between  the  determinant  and  the  permanent  (with  H. 
Mine),  lllinios  J.  Math.  5  (1961),  376-381. 
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